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TASK  OBJECTIVES 


There  are  two  broad  objectives  in  this  research  effort.  The  first  of  these 
involves  the  analysis  of  oceanic  data  sets  for  the  purpose  of  evaluating 
seafloor  noise  levels  vis-a-vis  stations  which  are  located  on  islands  and 
confirming  a  propagation  model  for  high  frequency  regional  phases  which 
was  developed  using  data  collected  during  the  DARPA-sponsored  NGENDEI 
Expedition.  The  second  objective  is  to  use  the  experience  with  regional  phases 
which  was  gained  in  studying  oceanic  Pn  and  Sn  in  analyzing  NORESS  array 
recordings  of  events  which  have  propagated  over  regional  distances  on  the 
continents. 

Specifically,  the  Work  Statement  includes  the  following: 

•  Obtaining  stable  estimates  of  the  short  period  (>0.1  Hz)  noise  at  island 
stations  (SNZO  and  GU MO)  for  comparison  with  typical  noise  levels  on 
the  seafloor. 

•  Examining  high  frequency  regional  phase  (Pn/Sn)  data  from  the 
Wake  Island  Array  for  comparison  with  the  Pn  propagation  model  and 
the  data  collected  in  the  southwest  Pacific  during  the  NGENDEI 
Expedition. 

•  Analyzing  data  from  the  NORESS  array  including  the  computation  of 
frequency-wavenumber  spectra  from  high  frequency  regional  phases 
such  as  Pg  and  Lg.  Comparing  the  results  with  computations  of 
synthetic  seismograms  and  frequency-wavenumber  spectra  in  order  to 
understand  more  fully  the  responsible  wave  propagation  phenomena. 
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TECHNICAL  RESULTS 


ABSTRACT 

In  this  report  we  present  the  results  of  several  investigations  dealing  with  two 

diverse  topics: 

•  We  compare  ambient  noise  levels  on  islands  with  those  present  at  and 

beneath  the  seafloor  to  determine  the  value  of  permanent  seafloor 
seismic  stations.  We  compare  signal  power  spectral  levels  expected  from 
earthquakes  with  ambient  noise  levels  at  these  locations  to  determine  a 

lower  magnitude  limit  of  detectable  events. 

•  We  study  the  nature  of  propagation  of  seismic  waves  through  the 

continental  and  oceanic  crusts.  We  investigate  the  influence  the 

regular  repetition  of  features  in  seismograms,  introduced  at  the  source 
and/or  during  propagation  by  layer  resonance,  has  on  the  spectrum  of 
the  recorded  coda. 


SUMMARY 

The  initial  goal  of  this  project  was  to  compare  ambient  noise  levels  on  islands 
with  those  present  at  and  beneath  the  seafloor.  The  motivation  for  this  work 
stems  from  concern  about  the  sparse  coverage  of  oceanic  areas  provided  by 
island  based  seismic  stations.  In  chapter  1  we  compare  noise  levels  recorded 
by  Seismic  Resarch  Observatory  (SRO)  stations  on  Guam,  the  north  island  of 
New  Zealand,  Taiwan  and  Easter  Island  and  conclude  that  in  the  frequency 
band  from  0.1  to  10  Hz  the  island  noise  levels  are  comparable  to  those  at  the 
seafloor,  recorded  by  Ocean  Bottom  Seismometers,  and  somewhat  higher  than 
those  recorded  beneath  the  seafloor  by  the  Marine  Seismic  System  (MSS).  We 
conclude  that  in  this  band  seafloor  stations  are  realistic  alternatives  when 
island  sites  are  not  available. 

In  chapter  2  we  extend  this  study  by  comparing  ambient  noise  levels  recorded 
on  and  beneath  the  seafloor  with  signal  levels  expected  at  an  epicentral  range 
of  30°  from  events  ranging  in  Mw  from  5.0  to  9.5.  We  find  that  at  this  range  in 

the  Pacific  events  with  Mw  as  low  as  5.5  should  be  detectable  by  seafloor 
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sensors.  Due  to  lower  noise  levels  in  the  Atlantic,  and  beneath  the  seafloor  in 
the  Pacific,  it  appears  that  the  threshold  for  sensors  deployed  in  these  regions 
is  as  low  as  5.0.  We  found  that  in  the  range  from  1  to  100  km  from  the 

shoreline  that  microseism  noise  levels  (in  the  band  near  250  mHz)  are 
inversely  correlated  with  distance  from  the  shoreline.  We  extend  this  effort, 
in  chapter  4,  with  a  comparative  analysis  of  ambient  borehole  noise  levels 
present  at  the  NORESS  array  in  Norway  and  in  Kazakhstan  U.S.S.R.  We  find 
that  the  noise  levels  at  the  NORESS  are  5  to  10  db  lower  than  those  in 
Kazakhstan  except  in  the  microseism  band  where  they  are  significantly 
higher.  We  conclude  this  is  probably  due  to  the  proximity  of  NORESS  to  the 
Ocean. 

In  chapters  3  through  6  we  consider  two  fundamentally  different  processes 
that  are  capable  of  impressing  a  regular  modulation  onto  the  frequency 
spectrum  of  recorded  coda.  The  first  process  we  consider  requires  the 
interaction  of  time  offset  wavefields  produced  by  sub-events  within  multiple- 
event  mine  explosions  (ripple-fired  quarry  blasts).  The  second  process 

involves  the  resonance  of  energy  in  low  velocity  strata. 

In  chapter  3  we  examine  recordings  ot  regional  single  and  multiple  event 

explosions  in  central  Asia.  We  consider  the  effect  that  the  practice  of  ripple 
firing,  or  the  detonation  of  several  explosions  closely  grouped  in  space  and 
time,  has  on  the  spectrum  of  the  resultant  seismic  coda.  We  predict 
theoretically,  and  observe  in  seismic  recordings,  that  this  practice  is  capable 
of  impressing  a  time-independent  modulation  on  the  spectrum  of  the  coda. 
Simple  events  do  not  seem  to  give  rise  to  this  time-independent  organization  of 
energy.  We  develop  an  algorithm  that  attempts  to  discriminate  between  mine 
explosions  and  single  event  explosions  by  seeking  time-independent  spectral 
features. 

In  chapter  4  we  extend  the  work  discussed  in  the  previous  chapter  by 

considering  data  collected  in  a  different  geographic  setting.  The  data  set 
considered  in  this  study  consists  of  recordings  of  earthquakes  and  mine 
explosions  made  by  the  NORESS  array  in  Norway.  Despite  low  digitization  rates 
at  NORESS  (40  samples/s)  the  two  types  of  events  are  easily  distinguished  using 
the  algorithm  developed  by  the  previous  study.  These  studies  show  that 
earthquakes  and  single  event  explosions  can  be  distinguished  from  mining 
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explosions  using  the  method  we  describe.  In  addition  the  algorithm  appears  to 
be  independent  of  geologic  setting  and  local  mining  practice. 

One  result  of  Scripps'  1983  S.W.  Pacific  Ngei.dei  expedition  was  the  Pn 
propagation  model  proposed  by  Sereno  and  Orcutt  (1985  a,b).  An  important 

aspect  of  this  model  is  a  coda  generation  mechanism  that  relies  solely  on  the 
resonance  of  energy  in  low  velocity  strata  and  does  not  invoke  scattering  by 
small  scale  inhomogeneities.  In  chapter  5  we  seek  to  reproduce  the  results  of 
Sereno  and  Orcutt  by  examining  high  frequency  recordings  of  earthquakes 
made  by  Hydrophones  in  the  Wake  Island  Hydrophone  array  and  by 
Hydrophones  and  Ocean  Bottom  Seismometers  deployed  during  the  Ngendei 
expedition.  Our  observations  are  essentially  consistent  with  those  made  by 
Sereno  and  Orcutt. 

In  chapter  6  we  report  on  research  carried  out  to  investigate  the  likely 

dominant  coda  generation  mechanism  present  in  the  continental  crust.  We 
use  two  datasets  for  this  purpose  -  recordings  of  mine  and  single-event 

explosions  made  in  central  Asia  as  well  as  recordings  of  mine  explosions  and 
earthquakes  made  by  the  NORESS  array  in  Norway.  This  work,  which 

primarily  involves  the  analysis  of  the  evolution  of  frequency  and 
wavenumber  spectra  as  a  function  of  time  in  seismic  coda,  is  currently  in 
progress  and  in  this  chapter  we  report  interim  results  and  plans  for  the 
continuation  of  this  research. 
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CHAPTER  1 


A  COMPARATIVE  STUDY  OF  ISLAND,  SEAFLOOR  AND  SUB¬ 
SEAFLOOR  AMBIENT  NOISE  LEVELS 

Michael  A.H.  Hedlin  and  John  A.  Orcutt 

Institute  of  Geophysics  and  Planetary  Physics  (A-025) 
Scripps  Institution  of  Oceanography,  La  Jolla  CA,  92093 

Bulletin  of  the  Seism.  Society  of  America  (1989)  vol.  79,  172-179 
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ABSTRACT 


A  study  of  seafloor  and  island  stations  shows  that  for  the  frequency  band  0.1  - 
10  Hz,  the  seismic  noise  levels  on  islands  are  comparable  to  the  levels  on  the 
seafloor.  The  microseism  peak  at  the  seafloor  appears  to  be  comparable  to  the 
highest  levels  observed  on  small  islands.  For  this  band,  seafloor  stations  are 
realistic  alternatives  when  island  sites  are  not  available. 

Seven  year  averages  of  the  ambient  noise  levels  recorded  by  Seismic  Research 
Observatory  (SRO)  stations  on  three  islands  (Guam  [GUMO],  Taiwan  [TATO]  and 
New  Zealand's  north  island  [SNZO])  are  compared  to  those  recorded  by  the 
International  Deployment  of  Accelerometers  (IDA)  station  on  Easter  Island  and 
on  and  beneath  the  ocean  floor  by  Ocean  Bottom  Seismometers  (OBS's)  and  the 
Marine  Seismic  System  (MSS)  deployed  in  a  south  Pacific  DSDP  drill  hole  at  23.8 
0  S.,  165.5  °  W  (Adair  et  al,1986).  From  0.3  to  2  Hz  the  SRO  displacement  power 
levels  fall  in  the  range  historically  observed  by  the  Scripps'  OBS's  (decreasing 
at  70  dB/decade  from  1  x  10^  nm^/Hz  at  0.3  Hz  to  1  nm^/Hz  at  2  Hz)  and  are  10  to 
15  dB  above  MSS  levels.  Above  2  Hz  it  appears  that  the  same  ratios  hold  (the 
SRO  power  levels  decrease  at  70  dB/decade  to  1  x  10'^  nm^/Hz  at  a  frequency  of 
10  Hz),  although  this  correlation  is  based  on  very  limited,  high  gain,  short 
period  data.  At  frequencies  below  0.3  Hz  the  SRO  noise  levels  peak  and  decrease 

to  approximately  2  x  10^  nm^/Hz  at  40  mHz.  The  noise  levels  recorded  at  Easter 

Island  are  somewhat  higher  (decreasing  at  70  dB/decade  from  1  x  10  ?  nm^/Hz 
at  0.2  Hz  to  1  nm^/Hz  at  10  Hz  and  to  1  x  10^  nm^/Hz  at  50  mHz).  At  the 

microseism  peak  near  0.2  Hz  the  MSS  levels  are  from  15  to  20  dB  higher  than 

observed  by  the  SRO  stations  and  equivalent  to  those  recorded  at  Easter  Island. 
There  appears  to  be  little  dependence  of  the  variance  in  noise  level  estimates 
on  frequency.  The  upper  95%  confidence  limit  generally  lies  10  dB  above  the 
average  noise  levels  for  all  island  stations. 

All  island  noise  level  curves  are  dominated  by  the  broad  double  frequency 
microseism  peak  centered  between  0.15  to  0.2  Hz.  The  single  frequency  peak 
ranges  from  absent  (Easier  Island)  to  discemable  (Guam  and  New  Zealand)  to 
obvious  (at  Taiwan),  The  center  frequency  of  this  peak  ranges  from  0.07  Hz  at 
Guam  and  New  Zealand  to  0.1  Hz  at  Taiwan.  We  speculate  that  the  increased 
amplitude  and  frequency  of  the  single  frequency  microseism  peak  is  due  to 
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the  interaction  between  the  shallow  continental  shelf  and  surface  gravity 
waves  and/or  the  presence  of  Taiwan  in  a  region  of  limited  fetch. 
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INTRODUCTION 


For  ihe  past  30  years,  seismologists  have  deployed  seismometers  on  and 
beneath  the  sea  floor.  The  primary  motivation  for  this  work  has  been  the 
desire  to  understand  ocean  crustal/lithospheric  structure  and  seismicity.  In 
addition,  however,  a  long  range  goal  of  this  research  has  been  to  use  these 
instruments  to  improve  the  global  coverage  of  seismometers.  To  demonstrate 
the  feasibility  of  long  term  ocean  bottom  deployments  i‘  must  be  shown  that 
noise  levels  at  the  seafloor  are  not  significantly  higher  than  those  commonly 
observed  on  land 

This  study  is  intended  to  assess  the  noise  level  change,  if  any,  that  can  be 
expected  by  moving  seismic  sensors  away  from  land.  Toward  this  end  this 
paper  compares  ambient  noise  levels  on  four  islands  with  those  at  and  below 
the  sea-floor.  In  addition,  the  island  noise  level  estimates  are  based  on  several 
noise  observations  spanning  a  number  of  years  so  that  an  estimate  of  the  long 
term  variance  of  these  levels  due  to  changing  meterological  and  cultural 
conditions  as  a  function  of  frequency  will  be  obtained.  The  locations  and 
relative  magnitudes  of  the  microseism  peaks  in  the  island  spectra  will  also  be 
discussed. 


THE  DATA 

The  SRO  data  used  in  this  study  were  taken  from  day  tapes  of  stations  deployed 
on  Guam,  Taiwan  and  the  north  island  of  New  Zealand,  (GUMO,  TATO  and  SNZO 
respectively,  see  table  1).  Short  and  long  period  SRO  data  (with  sampling  rates 
of  20  and  1  sample  per  second  respectively)  have  been  merged  to  provide 
power  estimates  over  the  range  from  40  mHz  to  10  Hz.  Only  the  vertical 
component  data  were  studied  since  no  horizontal  short  period  SRO  data  are 
collected  by  the  stations.  Each  power  estimate  is  based  on  averaging  10 
separate  samples  taken  at  random  from  4  to  7  year  time  spans  between  1976 
and  1983  so  that  a  wide  range  of  meteorological  conditions  were  incorporated. 
Prior  to  averaging,  each  trace  was  smoothed  so  that  the  final  variance  estimate 
would  be  largely  due  to  changing  meteorological  conditions  only. 
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Short  period  data  are  available  only  on  event  triggered  records  which  are 
relatively  infrequent.  Each  time  series  was  examined  on  a  visual  display 
before  transforming  to  ensure  that  data  corresponding  to  seismic  arrivals 
were  eliminated.  Each  long  and  short  period  sample  consisted  of  853  points 
padded  to  1024  points  with  zeros.  The  time  series  were  demeaned  and  then 
tapered  with  a  Hanning  window  prior  to  Fourier  transformation.  The  SRO 
seismometers  are  equipped  with  a  notch  filter  centered  at  roughly  0.17  Hz 
which  is  designed  to  reduce  the  effect  of  the  microseismic  noise  on  the  long 
period  channels.  This  notch  causes  a  sharp  peak  in  the  power  spectrum 
following  a  correction  for  the  instrument  response.  Since  this  peak  is 
unlikely  to  be  a  real  feature  of  the  data,  the  portion  of  the  power  spectra 
between  .15  and  .19  Hz  was  removed  and  the  plotted  levels  were  obtained  by 
interpolation.  The  power  levels  in  this  frequency  band  are  clearly  not  valid 
and  should  be  ignored.  For  a  discussion  of  SRO  instrumentation  the  reader  is 
referred  to  Peterson  et  al  (1976). 

The  IRIS/IDA  data  were  collected  by  the  vertical  component  of  the  Streckeisen 
seismometers  deployed  on  Easter  Island  (station  RPN,  see  table  1).  These  data 
were  treated  in  the  same  manner  as  the  SRO  data  although  only  9  intermediate 
period  (5  samples/s)  and  4  short  period  (20  samples/s)  samples,  taken  from 
June,  1987  to  March,  1988,  were  included  in  the  final  noise  level  estimate.  For 
a  discussion  of  the  instrument  refer  to  Wielandt  and  Streckeisen  (1982). 

The  MSS  data  were  collected  in  February  of  1983  during  Scripps'  Ngendei 
Expedition  and  the  noise  analysis  was  originally  published  by  Adair  et  al 
(1986).  The  MSS  was  deployed  at  23.8°  S,  165.5°  W  in  DSDP  drill  hole  595b  (leg 
91),  124  meters  below  the  sediment- water  interface  and  54  meters  into  the 
basement.  The  MSS  data  were  obtained  from  short  (40  samplcs/s)  and  mid¬ 
period  (4  samples/s)  channels  with  time  series  which  were  roughly  13  and  128 
seconds  long  respectively.  To  reduce  the  variance,  the  short  and  mid-period 
curves  are  based  on  the  average  of  10  spectra.  The  OBS  noise  level  bounds 
were  computed  from  25  noise  level  estimates  obtained  from  data  collected  at 
locations  and  dates  detailed  in  table  2.  Each  individual  estimate  was  made  after 
stacking  from  3  to  125  noise  level  curves  although  the  individual  time  series 
collected  were  too  short  to  obtain  high  resolution  at  frequencies  approaching 
0.2  Hz.  Both  MSS  and  OBS  power  estimates  have  been  corrected  for  instrument 
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response.  The  MSS  and  OBS  instruments  are  discussed  in  detail  in  Adair  et  al 
(1986)  and  Orcutt  et  al  (1987). 


COMPARISON  OF  NOISE  LEVELS 

The  four  island  stations  exhibit  similar  time-averaged  noise  structure  (figures 
1,2,3  and  4).  All  records  are  dominated  by  the  double  frequency  microseism 
peak  centered  between  0.15  and  0.2  Hz.  The  single  frequency  peak  ranges 
from  absent  (Easter  Island)  to  discemable  (Guam  and  New  Zealand)  to  obvious 
(at  Taiwan).  The  center  frequency  of  this  peak  ranges  from  0.07  Hz  at  Guam 
and  New  Zealand  to  0.1  Hz  at  Taiwan.  At  frequencies  above  the  main 
microseism  peak,  all  four  spectra  fall  off  at  roughly  70  dB/decad e.  In  all 
spectra  the  rate  of  decay  decreases  at  frequencies  above  2  Hz  and  the  SRO 
station  power  levels  actually  increase  above  5  Hz.  It  is  felt  that  this  apparent 
"levelling  off"  of  the  power  spectra  is  not  a  real  feature  of  the  data  but  is  due  to 
digitization  and  instrumentation  noise.  The  three  SRO  stations  used  in  this 
report  were  selected  because  they  are  located  on  islands  and  thus  are  close  to 
the  shoreline.  To  combat  the  resultant  high  noise  levels,  the  SRO  short  period 
station  gains  were  reduced  to  a  level  60  dB  lower  than  those  of  all  other  SRO 
stations,  (Peterson  et  al,1980).  As  a  result,  these  records  are  more  susceptible 
to  digitization  and  instrumentation  noise  at  high  frequencies.  In  the  early 
days  of  1976  at  GUMO,  however,  the  short  period  data  were  collected  with  high 
gain.  The  spectrum  from  12  January  1976  is  shown  in  figure  5.  At  frequencies 
above  2  Hz  the  power  spectral  levels  continue  to  decay  at  roughly  70 
dB/decade.  This  spectrum  is  clearly  more  representative  of  the  true  power 
level  in  this  frequency  band.  Because  the  average  long  term  power  levels  at 
all  stations  were  desired  it  was  obviously  not  possible  to  base  the  average 
power  level  estimates  solely  on  these  high  gain  data. 

At  frequencies  above  0.2  to  0.3  Hz  the  island  power  levels  fall  within  the  OBS 
noise  level  bounds  and  are  10  to  15  dB  above  MSS  levels.  At  frequencies  below 
0.2  Hz  the  MSS  levels  are  comparable  to  those  recorded  at  Easter  Island  and  10 
to  15  dB  above  noise  levels  recorded  on  the  3  SRO  islands.  This  enhanced  island 
noise  level  is  probably  a  reflection  of  the  oceanic  origin  of  the  microseisms. 
The  upper  95%  confidence  limit,  lying  two  standard  deviations  above  the 
average,  is  plotted  on  figures  1  to  4  with  the  average  power  level  curve.  The 
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lower  confidence  limit  was  excluded  since  it  was  negative  in  places.  These 
curves  track  each  other  fairly  well  and  indicate  that  each  frequency  band  has 
approximately  the  same  variance. 


THE  ISLAND  MICROSEISM  PEAKS 

The  traditional  explanation  of  the  origin  of  the  double  frequency  microseism 
peak  involves  the  non-linear  interference  of  gravity  waves  traveling  in 
opposite  directions,  Longuet-Higgins  (1950).  In  deep  water  this  non-linear 
interaction  is  effective  in  producing  acoustic  energy  at  twice  the  frequency  of 
the  interacting  gravity  waves.  This  acoustic  energy  subsequently  excites 
microseismic  activity  in  the  basement  rock  at  the  same  frequency.  The 
explanation  of  the  single  frequency  microseism  peak  is  fundamentally 
different.  Early  observations  of  this  peak  were  made  by  Wiechert  (1904), 

Oliver  and  Ewing  (1957)  and  Haubrich  et  al  (1963).  The  generally  accepted 
explanation  holds  that  when  surface  gravity  waves  propagate  through  shallow 
water  the  decaying  pressure  oscillations  they  cause  are  "felt"  directly  by  the 
substrate.  These  pressure  oscillations  thus  excite  seismic  energy  at  the  same 
frequency  as  the  source  surface  gravity  waves.  A  generally  accepted  rule  of 
thumb  holds  that  this  interaction  is  significant  when  the  water  depth  is  less 
than  60%  of  the  wavelength  of  the  surface  gravity  waves.  The  pressure 

fluctuations  caused  by  such  a  wave  on  the  bottom  is  greater  than  4.6%  of  the 

static  water  pressure.  At  0.07  Hz  a  typical  wavelength  is  on  the  order  of  350  m. 

This  crude  calculation  suggests  that  in  water  shallower  than  about  200  m, 
single  frequency  microseisms  can  be  excited. 

Adair  (1985)  discussed  previous  continental  microseismic  studies  which  have 
generally  detected  the  small,  single  frequency  microseism  peak  at  0.07  Hz  and 
the  larger  amplitude,  double  frequency,  peak  at  0.14  Hz.  The  single  frequency 
microseism  peaks  in  the  GUMO  and  SNZO  spectra  agree  with  each  other  and 
conform  to  these  continental  studies.  However,  at  TATO  the  fundamental 
microseism  peak  is  more  apparent  and  both  peaks  are  shifted  to  higher 
frequencies  (0.1  and  0.2  Hz).  The  single  frequency  peak  is  absent  from  the 
Easter  Island  data. 
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Plausible  explanations  for  these  differences  may  be  found  by  considering  the 
bathymetric  and  geographic  settings  of  the  four  islands.  Taiwan  is  located  on 
the  continental  shelf  of  China  and  is  ISO  km  from  the  mainland.  There  is  a 
broad  expanse  of  water  less  than  200  m  deep  in  the  vicinity  of  this  island  (the 
shallowest  water,  generally  less  than  100  m  deep,  lies  to  the  north  and  west). 
The  north  island  of  New  Zealand,  Guam  and  Easter  Island,  however,  are  located 
far  from  major  land  masses.  The  local  water  depths  increase  rapidly  away 
from  the  shore  to  greater  than  200  m  so  that  the  generative  area  for  primary 
microseisms  is  small  (relative  to  Taiwan).  This  is  especially  true  at  Guam  and 
Easter  Island.  It  seems  likely  that  the  increased  amplitude  of  the  primary 
microseismic  peak  at  TATO  relative  to  GUMO,  SNZO  and  especially  RPN  is  due  to 
the  broad  expanse  of  shallow  water  surrounding  Taiwan  which  facilitates  the 

direct  conversion  of  gravity  wave  pressure  fluctuations  to  seismic  energy  in 
the  basement. 

Hasselmann  and  Collins  (1967)  discussed  the  influence  of  shallow  water  on  the 
frequency  content  of  surface  gravity  waves.  As  these  waves  travel  in  shallow 
water,  the  long  period  components  that  interact  directly  with  the  sea  floor  are 
preferentially  dissipated.  In  the  vicinity  of  Taiwan  where  the  surface  gravity 

waves  most  likely  have  traveled  a  great  distance  over  the  shallow  basement 
before  approaching  the  island  it  seems  that  this  mechanism  would  strip  away 
the  low  frequencies  preferentially  and  cause  the  apparent  shift  of  the  single 

frequency  microseism  peak  to  a  slightly  higher  frequency. 

It  is  widely  known  that  the  frequency  (fm )  of  the  main  energy  peak  in  the 

spectrum  of  surface  gravity  waves  depends  on  the  fetch  and  the  local  wind 

speed  (commonly  represented  by  Ujg,  the  velocity  of  wind  10  meters  above  the 
sea  surface).  Hasselmann  et  al  (1973)  found  that  was  proportional  to  (U|q)' 
•34  in  the  North  Sea,  a  region  of  limited  fetch.  In  theory  fm  is  proportional  to 
(Uio)-10  in  unlimited  fetch.  The  waves  grow  until  the  phase  velocity,  which 
is  inversely  proportional  to  the  frequency,  equals  that  of  the  wind.  The 

proximity  of  Taiwan  to  mainland  China  restricts  the  local  wind  fetch  and  thus 

could  plausibly  explain  the  increased  frequency  of  the  single  frequency 
microseism  peak.  The  dependence  of  fm  on  the  wind  speed  itself  could  provide 

an  alternative  explanation  for  the  shift  in  Taiwan's  single  frequency 
microseism  peak  provided  that  the  long  term  average  wind  speed  was 
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significantly  lower  in  the  vicinity  of  Taiwan  than  near  the  other  islands. 
However,  wind  speed  data  published  by  the  U.S.  Naval  Weather  Service  indicate 
little  difference  between  the  average  scalar  wind  speeds  at  these  locations.  At 
Taiwan,  Guam,  New  Zealand  and  Easter  Island  the  average  wind  speeds  are  14.0, 
12.7,  14.5  and  13.0  knots. 


CONCLUSION 

We  have  computed  and  compared  estimates  of  noise  levels  on  four  islands,  on 
the  seafloor  and  beneath  the  seafloor.  Over  a  limited  frequency  band  the 
ocean  bottom  and  island  power  levels  are  comparable  and  thus  no  degradation 
of  the  signal  to  noise  ratio  should  be  expected  by  moving  sensors  from  land  to 
the  seafloor.  We  have  shown  that  over  this  frequency  band  a  sensor  located 
beneath  the  ocean  bottom  will  display  significantly  lower  noise  levels  than 
seismometers  deployed  on  the  ocean  bottom  or  on  islands.  This  is  true  except  at 
the  lowest  frequencies  (below  0.2  Hz).  These  observations  are  particularly 
encouraging  with  regard  to  sea  floor  stations  over  the  bulk  of  the  Pacific  and 
Indian  oceans  where  few  islands  are  available  for  station  siting. 

It  seems  plausible  that  the  increased  amplitude  and  frequency  of  the  single 
frequency  microseismic  peak  in  the  TATO  record  relative  to  GUMO,  SNZO  and 
RPN  is  due  to  the  interaction  between  the  shallow  continental  shelf 
surrounding  Taiwan  and  surface  gravity  waves  and/or  the  presence  of  Taiwan 
in  a  region  of  limited  fetch. 
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TABLE  AND  FIGURE  CAPTIONS 

Table  1.  The  geographic  coordinates  of  the  four  island-based 
seismometers  incorporated  in  this  study.  The  areas  of  the  islands  on  which 
these  seismometers  are  deployed  are  included. 

Table  2.  The  geographic  coordinates  and  dates  at  which  the  OBS  noise 
samples  used  in  this  study  were  collected. 

Figure  1.  Noise  power  spectral  density  in  displacement  units  for  the  SRO 

station  on  Guam.  The  lower  solid  curve  is  the  average  noise  level  obtained 
from  seven  years  of  data  while  the  upper  curve  is  a  95%  confidence  limit  on 
the  power.  The  frequently  dashed  curve  is  the  power  spectral  density  from 

the  Marine  Seismic  System  deployed  during  the  southwest  Pacific  Ngendei 
Expedition.  The  broadly  dashed  bounds  represent  upper  and  lower  limits  of 
noise  levels  observed  by  the  Scripps'  OBSs  over  a  twelve  year  period  in  the 
Pacific  Basin. 

Figure  2.  Noise  power  spectral  density  in  displacement  units  for  the  SRO 

station  on  the  north  island  of  New  Zealand.  The  other  curves  are  explained  in 
the  caption  for  Figure  1. 

Figure  3.  Noise  power  spectral  density  in  displacement  units  for  the  SRO 

station  on  Taiwan.  The  other  curves  are  explained  in  the  caption  for  Figure  1. 

Figure  4.  Noise  power  spectral  density  in  displacement  units  for  the 
IRIS/IDA  station  on  Easter  Island.  The  other  curves  are  explained  in  the 
caption  for  Figure  1. 

Figure  5.  Power  spectral  levels  in  displacement  units  for  GUMO  for  records 
in  1976  when  the  SRO  station  was  operating  at  higher  gain.  The  other  curves 
are  explained  in  the  caption  for  Figure  1.  Note  that  the  noise  levels  continue 
to  decay  at  higher  frequencies  indicating  that  the  higher  noise  levels 
observed  at  this  and  the  other  island  stations  are  probably  associated  with 
digitization  noise. 
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Table  1  Island  Data 


SITE 

LATITUDE 

LONGITUDE 

ISLAND  AREA 
(Square  Km) 

QUMO 

13°  35'  17"  N 

144°  51'  58"  W 

541 

SNZO 

41°  18'  36"  S 

174°  42’  16"  W 

115,000 

TATO 

24°  58’  32"  N 

121°  29’  18"  W 

36,000 

RPN 

27°  08'  37"  S 

109°  26'  10"  W 

118 
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Table  2  OBS  data 


NUMBER 

OF  SAMPLES 

LATITUDE 

LONGITUDE 

6 

41°  30'  N 

127°  20'  W 

2 

16°  O'  N 

145°  O’  W 

1 

29°  30’  N 

122°  O'  W 

3 

18°  O'  N 

145°  20'  W 

8 

16°  30*  N 

100°  30'  W 

4 

20°  50'  N 

109°  6'  W 

1 

23°  48*  S 

165°  30’  W 

DATE 

JUN  1976 
OCT  1977 
OCT  1977 
JAN  1976 
JUN  1977 
APR  1979 
FEB  1983 
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ABSTRACT 


Ambient  noise  levels  recorded  on  and  beneath  the  seafloor  are  compared  with 
the  signal  levels  expected  at  an  epicentral  range  of  30°  from  events  ranging 
in  Mw  from  5.0  to  9.5.  The  results  indicate  that  at  this  range  events  with  Mw  as 

low  as  5.5  should  be  detectable  by  sensors  on  the  seafloor  in  the  Pacific. 
Because  of  lower  ambient  noise  levels  present  in  the  Atlantic  and  beneath  the 
seafloor  in  the  Pacific,  events  with  Mw  as  low  as  5.0  should  be  detectable  in 

these  locations. 

Noise  level  data  collected  by  three  land-based  IDA/IRIS  sensors  (RPN,  ESK  and 
PFO)  ranging  in  distance  from  1  to  100  km  from  the  shoreline  are  compared. 
Low  frequency  noise  levels  (below  50  mHz)  are  not  sensitive  to  this  distance 
although  at  frequencies  above  50  mHz  noise  levels  are  clearly  inversely 
proportional  to  the  distance  from  the  shoreline.  At  250  mHz  noise  power 
spectral  levels  drop  40  dB  as  the  observation  point  varies  from  1  to  100  km 
from  the  coast. 
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INTRODUCTION 


It  is  obvious  that  the  global  coverage  which  land  based  seismometers  provide 
is  limited  since  71%  of  the  Earth's  surface  is  covered  by  water.  A  much  more 
uniform  coverage  is  feasible  if  instruments  can  be  deployed  on  the  seafloor  in 
environments  with  signal-to-noise  ratios  comparable  to  those  available  on 
land.  Many  technical  deployment  problems  have  been  overcome  during  the 
past  30  years  as  seismologists  have  deployed  seismometers  on  and  beneath  the 
sea  floor.  A  previous  paper  (Hedlin  and  Orcutt,  1989)  has  already  dealt  with  the 
question  of  the  relative  noise  levels  on  land  and  the  seafloor.  It  was  concluded 
that  over  the  frequency  band  from  40  mHz  to  10  Hz  the  ocean  bottom  and  island 
noise  levels  are  comparable.  In  this  paper  we  look  more  deeply  into  the 
acceptability  of  seafloor  noise  levels  by  comparing  them  with  signal  levels 
expected  from  earthquakes  of  varying  sizes  at  an  epicentral  range  of  30°. 

Although  seafloor  sites  could  potentially  provide  a  uniform  global  coverage, 
island  seismic  stations  still  play  an  important  role  in  the  global  seismographic 
network  since  they  provide  a  means  for  deploying  sensors  in  oceans  relatively 
inexpensively  and  with  relative  ease.  In  this  paper  we  will  also  examine  the 
dependence  of  island  and  continental  noise  levels  on  the  distance  from  the 
shoreline. 


THE  DATA 

The  sub-seafloor  noise  level  data  were  collected  by  the  Marine  Seismic  System 
(MSS)  in  February  of  1983  during  the  Scripps'  Ngendei  Expedition  and  the 
noise  analysis  was  originally  published  by  Adair  et  al  (1986).  The  MSS  was 
deployed  at  23.8°  S,  165.5°  W  in  DSDP  drill  hole  595b  (leg  91),  124  meters  below 
the  sediment-water  interface  and  54  meters  into  the  basement.  The  MSS  data 

were  obtained  from  short  (40  samples/s)  and  mid-period  (4  samples/s) 
channels  with  time  series  which  were  roughly  13  and  128  seconds  long 
respectively.  To  reduce  the  variance  of  the  power  spectral  estimates,  the  short 
and  mid-period  curves  are  based  on  the  average  of  10  spectra.  The  Ocean 
Bottom  Seismometer  (OBS)  noise  level  bounds  were  computed  from  25  noise 
level  estimates  obtained  from  data  collected  at  the  locations  and  dates  detailed 
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in  table  1.  Each  individual  estimate  was  made  after  stacking  from  3  to  125  noise 
level  curves  although  the  individual  time  series  collected  were  too  short  to 
obtain  high  resolution  at  frequencies  approaching  0.2  Hz.  Both  MSS  and  OBS 
power  estimates  have  been  corrected  for  instrument  response.  The  MSS  and 
OBS  instruments  are  discussed  in  detail  in  Adair  et  al  (1986)  and  Orcutt  et  al 
(1987). 

The  data  for  the  high-frequency  seafloor  pressure  spectra  were  collected  on 
two  separate  occasions.  The  Pacific  data  were  collected  in  May  1988  during  the 
Scripps'  Nachos  Expedition  350  km  west  of  San  Diego  in  3850  m  of  water.  A  low 
frequency  hydrophone  (Cox  et  al,  1984)  mounted  on  the  side  of  an  OBS  was 
sampling  at  32  samples/s  and  the  recorder  collected  15  minute  records.  The 
28,800  point  time  series  were  divided  into  30%  overlapping  8192  point  windows. 
Each  series  was  Hanning  tapered  prior  to  Fourier  transformation  and  then 
corrected  for  the  instrument  response.  The  final  spectral  estimate  was 
obtained  from  a  stack  of  the  resulting  spectra.  The  Atlantic  seafloor  pressure 
spectrum  was  computed  from  data  collected  by  the  Woods  Hole  Oceanographic 
Institution  during  the  HEBBLE  expedition  in  October  of  1985.  The  HEBBLE 
pressure  data  were  collected  by  a  long  period  transducer  (Cox  et  al,  1984) 
deployed  at  4820  m  depth  over  the  Nova  Scotia  Rise.  The  hydrophone  data  were 
sampled  at  one  sample/s  during  6  hr  time  periods  each  day.  To  compute  the 
power  spectral  estimate,  the  data  were  treated  in  a  similar  manner  as  the 
Atlantic  high  frequency  data,  the  only  difference  being  that  4n  prolate 
windows,  rather  than  Hanning  windows,  were  used  to  taper  each  sub¬ 
sequence. 

The  land-based  data  were  collected  exclusively  at  IDA/IRIS  stations.  Low  (1 
sample/s)  and  high  frequency  (20  samples/s)  spectral  estimates  have  been 
merged  to  provide  the  broadband  spectra  in  figure  3.  Each  spectrum  is  the 
result  of  stacking,  in  the  frequency  domain,  several  spectral  estimates 
calculated  from  distinct  portions  of  the  original  noise  sample.  The  very  low 
frequency  seafloor  pressure  data  were  obtained  from  figure  6  of  Filloux  (1980). 
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COMPARISON  OF  SEAFLOOR  SIGNAL  AND  NOISE  LEVELS 


Figure  1  illustrates  three  seafloor  pressure  power  spectra  spanning  the  band 
from  0.7  pHz  to  20  Hz.  The  solid  and  finely-dashed  curves  represent  data 
collected  in  the  Pacific  while  the  data  for  the  lower  power,  coarsely-dashed 
curve  were  collected  in  the  Atlantic.  Power  spectral  levels  in  the  lowest  band, 
from  0.7  pHz  to  20  mHz,  decrease  monotonically  at  25  dB/decade  except  in  the 
band  from  10  to  200  pHz  where  several  tidal  lines  are  recorded.  The  microseism 
band  in  both  spectra  is  dominated  by  the  double  frequency  microseism  peaks 
above  200  mHz.  Between  the  microseism  peaks  and  2  Hz,  power  spectral  levels 
decrease  at  40  dB/decade. 

In  figure  2  the  high-frequency  data  in  Figure  1  have  been  converted  to 
acceleration  units  and  superimposed  on  curves  which  illustrate  the 
acceleration  expected  at  an  epicentral  range  of  30°  from  earthquakes 
spanning  a  range  of  magnitudes.  OBS  and  MSS  data  have  been  plotted  as  well. 
The  conversion  from  pressure  to  amplitude  was  made  using  a  simple  acoustic 
plane  wave  relationship  assuming  vertical  incidence  of  seismic  energy  and 
continuity  of  vertical  traction  and  displacement.  Although  the  validity  of  this 
conversion  is  clearly  restricted,  we  expect  that  the  predicted  accelerations  are 
probably  realistic  estimates.  The  validity  of  the  noise  levels  can  only  be 
verified  through  experimentation.  This  figure  is  very  encouraging  and 
indicates  that  with  the  exception  of  the  noisiest  OBS  sites  all  events  with  Mw  of 

5.5  or  greater  lead  to  motions  above  the  expected  noise  levels  at  a  range  of  30°. 
Such  earthquakes  or  explosions  should  be  detectable  with  seafloor  sensors. 
Using  hydrophones  deployed  in  the  Atlantic  and  the  MSS  buried  in  the 
substrate  beneath  the  Pacific,  where  ambient  noise  levels  appear  to  be 

somewhat  lower,  earthquakes  with  Mw  as  low  as  5.0  should  be  detectable  at  a 

range  of  30°.  This  conclusion  is  consistent  with  the  results  of  Orcutt  and 

Jordan  (1985).  Using  data  collected  by  the  MSS  during  the  Ngendei  Expedition 
they  found  that  no  events  with  m^  less  than  5.1  were  detected  at  ranges 

greater  than  30°.  At  lesser  distances,  where  high  frequency  signals  propagate 
through  the  high-Q  lithospheric  waveguide,  they  found  that  the  detection 
threshold  was  far  lower.  Shallow-focus  events  at  distances  on  the  order  of  15° 

with  mjj  as  low  as  3.7  could  be  detected  by  the  MSS.  At  the  noisiest  sites, 
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depending  on  the  frequency  band  of  interest,  the  smallest  events  detectable 
range  in  size  from  Mw  of  6  to  8. 

THE  DEPENDENCE  OF  LAND-BASED  NOISE  LEVELS  ON  THE  DISTANCE 

FROM  THE  SHORELINE 

Figure  3  shows  power  spectra  of  noise  data  obtained  from  the  IRIS/IDA 
seismographic  stations  on  Easter  Island  (RPN),  one  km  from  the  coastline,  at 
Eskdalemuir  (ESK),  40  km  from  the  Irish  Sea  and  at  Pifion  Flat  (PFO),  some  100 
km  from  the  California  coast.  The  dashed  lines  show  the  level  of  minimum  and 
maximum  ambient  noise  at  quiet  continental  locations.  The  high  frequency 
portion  of  the  spectrum  for  the  RPN  station  was  derived  from  data  collected 
during  the  most  noisy  period  experienced  during  a  one-year  period;  it  should 
be  taken  as  representative  of  the  maximum  high  frequency  noise  at  an  island 
site.  The  spectra  demonstrate  that  for  periods  longer  than  the  microseisms, 
there  is  only  a  small  decrease  in  noise  as  distances  to  the  coastline  increase. 
Island  stations  should  provide  relatively  low  noise  data  in  this  frequency 
range.  At  frequencies  higher  than  the  microseisms,  there  is  considerable 
reduction  in  noise  as  the  distance  from  the  coastline  increases.  At  250  mHz  the 
noise  levels  drop  40  dB.  Factors  such  as  coastal  shoaling,  barrier  reef 
coverage,  and  surf  conditions  can  be  expected  to  cause  large  variations  in  the 
ambient  noise  in  this  frequency  range. 

CONCLUSION 

It  appears  that  most  locations  on  the  seafloor  are  sufficiently  quiet  to  permit 
the  detection  of  moderate-sized  earthquakes  and  explosions  at  teleseismic 
ranges.  This  result  is  particularly  encouraging  in  view  of  the  need  for  an 
increased  uniformity  of  the  global  coverage  of  seismic  sensors. 

The  acceptability  of  island  sensor  deployments  depends  on  the  frequency  band 
of  interest  and  on  the  distance  of  the  sensor  from  the  shoreline  (and  thus  the 
size  of  the  island).  It  appears  that  low  frequency  (below  50  mHz)  noise  levels 
are  insensitive  to  the  distance  of  the  sensor  from  the  shoreline.  However,  in 
high-frequency  studies  it  appears  that  sites  on  large  islands  relatively  far 
from  the  shoreline  will  have  superior  noise  levels. 
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TABLE  AND  FIGURE  CAPTIONS 

Table  1.  The  geographic  coordinates  and  dates  at  which  the  OBS  noise 

samples  used  in  this  study  were  collected. 

Figure  1.  Noise  power  spectral  density  in  pressure  units  for  three  seafloor 

sensors.  The  data  for  the  solid  and  finely-dashed  curves  were  collected  in  the 
Pacific  Ocean,  off  the  Coast  of  California  (in  3850  m  of  water)  and  off  the  Gulf 
of  California  (at  3210  m  depth)  respectively.  The  data  for  the  coarsely-dashed 
curve  were  collected  in  the  Atlantic  at  the  Nova  Scotia  Rise  in  5000  m  of  water. 

Figure  2.  Acceleration  amplitudes  expected  at  an  epicentral  range  of  30° 

from  earthquakes  ranging  in  Mw  from  5.0  to  9.5.  The  data  for  the  light-solid 
curve  were  collected  in  the  Pacific  Ocean,  off  the  coast  of  California  at  3850  m 
depth.  The  two  heavy-solid  curves  represent  the  extreme  noise  level  bounds 

detected  at  25  sites  in  the  Pacific  (see  table  1).  The  frequently  dashed  curve 
represents  the  noise  level  recorded  by  the  Marine  Seismic  System  deployed 
during  the  southwest  Pacific  Ngendei  Expedition.  The  data  for  the 
infrequently  dashed  curve  were  collected  in  the  Atlantic  at  the  Nova  Scotia 
Rise  in  5000  m  of  water.  The  amplitudes  plotted  are  1.25  times  the  RMS  values 
over  1/3  octave  bandwidths. 

Figure  3.  Noise  recordings  made  on  an  island  (RPN-Easter  Island)  compared 

with  measurements  made  on  continental  sites  located  near  the  ocean  (ESK- 
Eskdalemuir;  PFO-Pifion  Flat  (Observatory).  Note  that  the  high  frequency 
content  of  the  microseism  noise  decreased  as  the  distance  from  the  shore 
increases.  The  amplitude  of  the  microseism  peak  also  decreases  in  the  same 
sense  as  the  high  frequencies. 
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Table  1  OBS  data 


NUMBER 
OF  SAMPLES 


LATITUDE 

LONGITUDE 

41°  30’  N 

127°  20'  W 

16°  O'  N 

145°  O'  W 

29°  30'  N 

122°  O’  W 

CD 

o 

O 

z 

145°  20'  W 

16°  30'  N 

100°  30'  W 

20°  50'  N 

109°  6'  W 

23°  48’  S 

165°  30’  W 

DATE 

JUN  1976 
OCT  1977 
OCT  1977 
JAN  1976 
JUN  1977 
APR  1979 
FEB  1983 


27 


(zj-i | Dosod )  JaMod  ajnssaj 

28 


10- 


29 


oi-Ol 


Figure 


CHAPTER  3 


THE  TIME-FREQUENCY  CHARACTERISTICS  OF  QUARRY 
BLASTS  AND  CALIBRATION  EXPLOSIONS  RECORDED  IN 

KAZAKHSTAN  U.S.S.R. 

Michael  A.  H.  Hedlin 
J.  Bernard  Minster 
John  A.  Orcutt 

Institute  of  Geophysics  and  Planetary  Physics  (A-025), 
Scripps  Institution  of  Oceanography, 

University  of  California  at  San  Diego,  La  Jolla  CA  92093 

accepted  for  publication  in  the  Geophysical  Journal 


31 


SUMMARY 


In  this  paper  we  consider  two  fundamentally  different  processes  that  can  be 
responsible  for  the  organization  of  energy  in  seismic  coda  into  discrete  time- 
independent  frequency  bands.  One  process  involves  the  resonance  of  energy 
in  low  velocity  horizons  and  the  other  requires  the  interaction  of  time  offset 
wavefields  produced  by  sub-events  within  multiple-event  mine  explosions 
(ripple-fired  quarry  blasts). 

We  examined  data  collected  by  high  frequency  seismometers  in  Kazakhstan 

U.S.S.R.  and  observed  regular  time-independent  spectral  modulations  in  coda 
resulting  from  events  strongly  suspected  to  be  Soviet  quarry  blasts,  but  not  in 

the  coda  from  single  event  calibration  explosions  detonated  at  similar  ranges. 

We  conclude  these  modulations  are  a  source  effect  and  due  to  ripple-firing. 
This  modulation  is  independent  of  the  source-receiver  azimuth  and  we  infer 
that  the  spatial  array  of  sub-shots  in  each  event  must  be  small.  We 
demonstrate  that  simple  linear  superposition  theory  can  be  used  to  reproduce 
effectively  the  spectral  modulation  observed  in  real  quarry  blasts. 

On  the  basis  of  these  observations  we  attempt  to  discriminate  between  the  two 
types  of  events  using  a  spectral  pattern-based  algorithm  that  seeks  time- 
independent  features.  We  consider  the  detrimental  effect  that  resonant 

energy  in  low  velocity  horizons  can  have  on  the  successful  application  of  our 
algorithm. 
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INTRODUCTION 


In  theory  any  process  that  introduces  a  regular  repetition  into  a  seismic 
wavetrain  will  impress  upon  the  spectrum  of  that  energy  a  regularly  spaced 
modulation.  Superposed  regular  modulations  were  observed  in  spectra 
computed  from  Pn  and  Sn  coda  collected  in  the  S.W.  Pacific  in  1983  during 
Scripps'  Ngendei  expedition.  These  were  explained  by  Sereno  and  Orcutt  (1985 
a  and  b)  as  being  due  to  the  resonance  of  energy  in  the  water  and  sedimentary 
columns  near  the  receiver.  This  mechanism  has  provided  an  elegant 
alternative  explanation  to  the  scattering  theory  previously  held  by  many  as 
the  dominant  source  of  coda  generation  in  the  oceanic  lithosphere. 

We  have  recently  observed  very  similar  spectral  modulation  in  data  collected 
in  the  Soviet  republic  of  Kazakhstan  in  1987.  This  data  set  includes  calibration 
events,  where  the  source-time  function  is  presumably  relatively  simple,  and 
many  other  regional  events,  a  large  fraction  of  which  are  likely  mine 
explosions.  This  unusual  spectral  modulation  is  only  observed  in  the  latter 
events  and  thus  we  suspect  that  it  is  not  a  propagation  effect,  as  was  the  case  in 
the  Ngendei  data  set.  For  example,  Baumgardt  and  Ziegler  (1988)  discuss  a 
commonly  used  technique  in  mining  known  as  ripple-firing.  In  this  practice 
a  number  of  sub-explosions  are  staggered  in  time  and  space  to  reduce  ground 
motions  in  areas  proximal  to  the  mine  as  well  as  to  enhance  the  fracturing  of 
the  rock.  The  most  plausible  explanation  of  the  spectral  modulation  we 
observe  in  the  Soviet  data  is  that  it  is  a  source  effect  which  stems  from  the 
highly  repetitive  nature  of  ripple-fired  mine  explosion  source-time  functions. 

Although  the  physics  behind  these  two  types  of  modulations  is  fundamentally 
different,  the  effect  they  have  on  the  spectra  of  the  recorded  energy  is 
strikingly  similar.  We  predict  mathematically  that  both  processes  are  capable 
of  impressing  a  time-independent  modulation  on  the  spectra  of  seismic 
wavefields. 

We  feel  the  problem  of  discriminating  ripple-fired  mine  explosions  (quarry 
blasts)  from  underground  nuclear  explosions  has  become  increasingly 
important  and  deserves  close  scrutiny.  As  discussed  by  Stump  and  Reamer 
(1988),  a  reduced  Threshold  Test  Ban  Treaty  would  bring  the  magnitude  of  the 
largest  allowable  nuclear  explosions  down  to  that  of  large  ''engineering" 
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explosions  otherwise  known  as  quarry  blasts.  The  problem  of  discriminating 

quarry  blasts  from  other  events  is  not  a  new  one,  and  has  been  investigated  by 
several  other  authors.  Aviles  and  Lee  (1986)  considered  quarry  blasts  and 
earthquakes  on  the  west  coast  of  the  United  States  and  found  they  had  success 

in  discriminating  between  the  two  types  of  events  by  computing  ratios  of  the 

spectral  energy  at  1.5  and  12  Hz.  They  found  quarry  blasts  to  be  relatively 

deficient  in  high  frequency  energy  and  proposed  a  preferred  attenuation  of 

high  frequency  energy  near  the  quarry  blast  source  by  near-surface 

fractures.  Although  this  method  may  be  valid  in  the  western  U.S.  it  seems  that 
its  success  would  be  highly  dependent  on  mining  practice  and  the  geology  of 

the  region  in  which  the  discriminant  is  being  used.  As  we  discuss  later,  the 

ratios  of  power  in  any  two  frequencies  is  highly  dependent  on  the  geometry 

and  timing  of  the  sub-shots  in  the  quarry  blast  sequence. 

Many  authors  have  performed  cepstral  analyses  of  multiple  source  events. 

The  cepstrum,  which  is  the  Fourier  transform  of  the  log  of  the  spectrum 

(Tribolet,  1979),  is  sensitive  to  regular  spectral  modulation  in  the  coda  of 

multiple  source  events.  Baumgardt  and  Ziegler  (1988)  showed  that  cepstra  of 

ripple-fired  sources  display  power  highs  at  certain  times  dependent  on  the 
timing  of  the  shots  in  the  quarry  blast  sequence.  A  weakness  of  studies  which 
employ  the  cepstrum  for  discrimination  purposes  is  that  simple  events 

(earthquakes  and  single  shot  explosions)  can  possess  significant  cepstral 
structure.  The  cepstra  computed  by  Baumgardt  and  Ziegler  (1988)  from  the 
coda  of  some  earthquakes  illustrate  this  point.  The  problem  of  discrimination 

becomes  somewhat  arbitrary  when  it  must  be  decided  how  much  structure  is 
sufficient  to  identify  an  event  as  a  quarry  blast. 

Gupta  et  al  (1984)  attempted  to  discriminate  between  single  and  multiple-event 
explosions  and  earthquakes  and  explored  several  methods  which  were  based 
on  the  expected  differences  between  the  distribution  of  energy  in  Pg  and  Lg 

phases  at  a  variety  of  frequencies.  The  expected  differences  they  discussed 
were  mainly  due  to  the  proximity  of  explosions  to  the  free  surface  and  their 
omnidirectionality  (resulting  in  their  relative  inability  to  generate  SH 
energy).  They  found  that  each  of  the  proposed  discriminants  taken  separately 
had  a  high  failure  rate  and  thus  it  was  necessary  to  combine  all  of  them  into  a 
multivariate  discriminant.  They  observed  large  differences  in  the  Pg  and  Lg 
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phases  excited  by  two  proximal  nuclear  explosions  and  attributed  the  wide 
variability  in  their  results  to  the  sensitivity  of  the  Pg  and  Lg  phases  to  near 
source  effects. 

We  feel  that  a  successful  discriminant  must  be  relatively  independent  of  the 
variability  of  seismic  phases,  of  mining  practice  and  geologic  setting.  In  this 
paper  we  suggest  that  the  property  that  best  distinguishes  multiple-shot 
explosions  from  all  other  seismic  events  is  the  persistence  through  time  of 
prominent  spectral  features  in  the  onsets  and  coda  produced  by  these  complex 
events.  We  propose  a  time  versus  frequency  pattern-based  discriminant 
which  seeks  long-lived  spectral  features  in  seismic  coda.  The  concept  of 
exploiting  the  time-independence  of  the  spectral  modulation  induced  by 
ripple-firing  the  source  is  not  entirely  new  and  was  previously  suggested  by 
Bell  (1977). 


SPECTRAL  MODULATIONS  IN  SEISMOGRAMS 

Sereno  and  Orcutt  (1985b)  demonstrate  the  significant  influence  of  low 
velocity  horizons  in  the  oceanic  environment  on  coda  duration  and  amplitude. 
This  layer  resonance  enriches  the  spectrum  of  the  coda  at  certain  equispaced 
eigenfrequencies,  determined  by  the  layer  thickness  and  velocity. 
Considering  a  relatively  simple  case  in  which  a  wavelet  w(t)  is  reverberating 
in  a  single  low  velocity  surface  layer,  the  linear  sum  of  all  reverberations  x(t) 
equals: 


x(t)  =  w(t)  * 


•  e  "’"  •  H(t) 


-iret/T 


(1) 


The  time  between  successive  reverberations  is  represented  by  T  and  to  a  good 
approximation  equals  the  two-way  vertical  travel  time  in  the  horizon.  The 
"shah"  function  is  tapered  exponentially  to  simulate  the  effect  of  attenuation 
and  is  multiplied  by  the  Heaviside  step  function  to  exclude  any  acausal  energy. 
The  final  exponential  is  required  by  the  n  phase  shift  incurred  by  reflection  at 
the  free  surface. 


35 


By  Fourier  transforming  the  preceding  equation  the  authors  found  that  the 
spectrum  of  the  original  wavelet  is  multiplied  by  an  infinite  set  of  staggered 
tapering  functions  which  decay  at  a  rate  controlled  by  and  proportional  to  the 
attenuation  parameter  a.  The  amplitude  spectrum  of  the  sum  of  wavelets  is 
given  by: 


1 

a  -  i2tif 

IX(f)  1  =  IW(f)J  • 

III(fT  +  f)  * 

(a2  -t-  47t2f2) 

(2) 


The  thicker  -  or  equivalently,  the  slower  -  the  horizon  the  more  closely  spaced 
the  spectral  highs  will  be.  In  the  simple  case  where  one  horizon  is 
responsible  for  the  bulk  of  the  resonance  and  the  receiver  is  a  significant 

distance  from  the  source  (so  that  the  incident  energy  has  high  phase  velocity) 

the  period  of  the  modulation  observed  on  vertical  or  horizontal  component 
sensors  is  largely  independent  of  time  in  the  coda  (Hedlin,  Orcutt  and  Minster, 
1988).  Sereno  and  Orcutt  (1985b)  point  out  that  in  situations  where  two  or 
more  different  modulations  are  superposed,  either  due  to  more  than  one 

resonating  horizon  or  the  combination  of  compressional  and  shear  resonance 
in  the  same  horizon,  the  interference  of  the  competing  modulations  can 
produce  time-dependent  behavior,  most  likely  manifested  as  an  evolution  to 

lower  frequency  content  as  time  progresses. 


We  now  consider  the  theory  underlying  a  distinctly  different  origin  of 
spectral  modulation,  one  originating  at  the  source  rather  than  through 
propagation.  We  consider  a  simple  model  of  a  ripple-fired  mine  explosion 
(quarry  blast)  in  which  all  sub-shots  occur  at  the  same  point  in  space  and  are 
offset  from  each  other  by  equal  time  spans  T.  In  addition  we  assume  that  all 
wavelets  produced  by  each  sub-shot  w(t)  are  identical  and  superpose  linearly. 
Approaching  the  problem  in  the  manner  of  Sereno  and  Orcutt  (1985b)  we  can 
compute  the  seismogram  x(t)  produced  by  the  entire  suite  of  sub-shots  by 
convolving  the  wavelet  (coda)  produced  by  each  identical  sub-explosion  by 
the  finite  duration  "shah"  function: 


x(t)  =  w(t)  * 


><T> '  B%> 


(3) 
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The  shah  function  is  multiplied  by  the  boxcar  function  B(t)  which  lasts  D 
seconds,  the  duration  of  the  entire  suite  of  sub-shots  in  the  quarry  blast.  By 

Fourier  transforming  (3)  we  find  that  the  spectrum  of  the  entire  seismogram 

equals  that  of  an  individual  shot  multiplied  by  an  infinite  set  of  equispaced 

"sine"  functions. 


X(f)  =  W(f)  • 


III(fT)  *  sin^ffQ) 


(4) 


Although  the  physics  behind  the  two  processes  that  lead  to  modulations  are 
fundamentally  different,  the  mathematics  that  describe  them  at  this  simple 
level  are  extremely  similar.  The  principal  difference  between  the  two  models 
arises  from  the  shah  function  which  convolves  the  starting  wavelet  and  is 
truncated  by  the  boxcar  describing  the  abrupt  beginning  and  end  of  the 
quarry  blast  source.  In  the  Pn  propagation  problem  the  overall  modulation 

arises  from  a  smoother  function  which  reflects  the  loss  of  energy  due  to 
attenuation.  Figure  1  displays  the  modulation  caused  by  linearly  superposing 

the  wavelets  produced  by  1  to  5  sub-shots  located  at  the  same  point  in  space  and 
delayed  100  ms  from  each  other.  The  amplification  at  the  preferred 
frequencies  equals  the  number  of  sub-shots.  The  degree  to  which  the  energy 
is  confined  to  the  preferred  frequencies  is  directly  dependent  on  the  duration 
of  the  quarry  blast  set.  By  inspection  of  equation  (4)  one  can  see  that  the 
smaller  the  spacing  between  successive  shots,  the  more  broadly  spaced  the 

modulation  in  the  frequency  domain  will  be.  The  frequency  spacing  is 
equivalent  to  the  inverse  of  the  shot  spacing  and,  as  is  displayed  in  figure  1,  a 
100  ms  spacing  results  in  a  frequency  spacing  of  10  Hz.  In  this  idealized  case 
the  modulation  is  independent  of  time. 

More  general  mine  explosions,  involving  arbitrary  offsets  in  time,  can  be 

simulated  in  a  similar  manner.  The  theory  has  been  dealt  with  by  Baumgardt 
and  Ziegler  (1988)  but  is  similar  to  that  described  above  provided  the  scatter  of 
offset  times  is  not  too  great.  Again  linear  superposition  is  assumed. 
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x(t)  = 


w(t-T.) 


(5) 


The  absolute  time  of  each  of  the  n  sub-events  (relative  to  the  first  event) 
equals  Tj.  After  Fourier  transformation  we  find  that  the  amplitude  spectrum  of 
the  sum  of  wavelets  equals: 


X(f)  =  W(f)  • 


cos(2rcfT.) 


i=l 


V 


sin(27tfT.) 


J 


(6) 


The  energy  from  one  sub-event  is  multiplied  by  a  sum  of  sines  and  cosines, 
and  is  still  organized  preferentially  into  time-independent  frequency  bands 

provided  that  the  scatter  of  time  delays  from  a  constant  value  is  not  too  great. 

Figure  2  displays  the  modulation  pattern  expected  from  an  event  consisting  of 
50  sub-shots.  The  time  delays  between  adjacent  events  are  distributed 
normally  with  a  mean  value  of  63  ms  and  a  variance  of  6.3  ms.  The  power  is 
concentrated  near  the  simply  predicted  frequencies  of  16,  32  and  48  Hz, 
although  as  the  frequency  of  the  overtone  increases  so  does  the  scatter.  This 
model  is  undoubtedly  more  realistic  than  the  first.  Although  the  mining 
engineers  may  want  their  sub-shots  to  adhere  to  a  regular  time-spacing  there 
may  be  considerable  deviation  from  this  in  practice  (Stump  and  Reamer,  1988). 
Using  high  speed  photographic  observations  of  the  sub-shot  blast  times  in 
quarries  in  the  eastern  U.S.  the  authors  observed  that  deviations  between  the 
intended  and  actual  times  were  as  high  as  34%.  In  addition,  no  true  quarry 
blast  consists  of  a  number  of  discrete  events  all  located  at  the  same  point  in 

space.  True  quarries  employ  a  spatial  array  of  sub-shots  in  blasting.  These 

spatial  offsets  are  equivalent  to  apparent  offsets  in  time  that  depend  on  the 
slowness  of  the  energy  being  considered  and  the  azimuth  from  the  mine  to  the 
receiver.  They  will,  along  with  the  true  scatter,  introduce  a  deviation  from  a 
regular  time  spacing  but  will  not  altogether  destroy  the  organization  of  the 
energy  into  time-independent  bands  provided  that  the  shot  array  dimensions 
are  not  too  great.  Although  this  is  quite  easily  dealt  with  (Smith,  1989),  as  will 
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be  discussed  later  it  appears  that  we  do  not  need  this  refinement  to  interpret 
the  events  we  examine  in  this  study. 


THE  DATA  SET 

The  data  used  for  this  study  were  collected  in  1987  in  the  Soviet  republic  of 
Kazakhstan  roughly  200  km  from  the  Semipalatinsk  underground  nuclear  test 

site.  Three  receiver  sites  were  established  in  the  spring  of  1987  in  a 
triangular  array  near  the  towns  of  Bayanaul,  Karkaralinsk  and  Karasu  (see 

figure  3).  The  local  crust  varies  in  thickness  from  roughly  40  km  (near 
Karasu)  to  SO  km  near  Karkaralinsk  and  Bayanaul  (Balyaevsky  et  al,  1973; 

Leith,  1987).  All  three  sites  were  installed  on  granitic  intrusions  of  Permian  to 
Triassic  age  and  consisted  of  surface  and  borehole  instruments  (deployed  at  99, 
66  and  101  meters  respectively).  Although  establishing  the  receivers  on 

granitic  outcrops  should  have  eliminated  the  problem  of  near  receiver 
reverberations,  the  site  at  Karasu  was  clearly  contaminated  by  an  anomalous 
site  response  (Eissler  et  al,  1988)  and  the  data  from  this  site  have  been  excluded 
from  this  study. 

Teledyne  Geotech  54100  seismometers  sampling  at  250  s'1  with  a  flat  response 
to  velocity  between  0.2  and  100  Hz  were  deployed  in  each  of  the  boreholes 

(Berger  et  al,  1987)  .  The  surface  data  used  in  this  study  were  collected  by  GS- 
13  seismometers  deployed  in  shallow  vaults,  which  were  also  sampled  at  250  s* 
1.  These  instruments  have  a  flat  response  to  ground  velocity  between  1  and  80 
Hz.  All  seismometers  recorded  data  only  after  triggering. 

Three  calibration  events  detonated  in  September  of  1987  provide  the  only  well 

constrained  events  in  this  data  set.  They  varied  in  size  from  10  (events  1  and 

3)  to  20  tons  (event  2)  and  ranged  from  157  to  254  km  from  the  receivers  (see 

figure  1).  Unlike  events  1  and  3,  event  2  was  not  fully  contained  underground. 
These  events  were  estimated  by  Eissler  et  al  (1988)  to  have  equivalent  seismic 
moments  on  the  order  of  10^2  to  10^  Nm.  Since  event  3  was  detonated  at  the 
same  site  as  event  1  and  occurred  simultaneously  with  the  arrival  of  a 
teleseism  originating  in  the  MacQuarie  Islands,  it  has  not  been  used  in  this 
study. 
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The  region  surrounding  the  seismic  network  experiences  a  low  level  of 

seismicity  (Leith,  1987)  and  is  considered  to  be  tectonically  stable.  We  believe 
that  the  bulk  of  the  events  recorded  during  the  operation  of  the  network  are 

man-made  and  are  likely  mine  explosions.  The  events  considered  in  this  paper 
range  from  105  to  264  km  from  the  receivers.  Compelling  support  for  the 

identification  of  one  of  these  events  (event  i,  figure  3;  table  1)  as  a  quarry 

blast  comes  from  a  French  SPOT  photo  provided  by  Dr.  Cliff  Thurber  at  the  State 

University  of  New  York,  Stony  Brook.  The  photo  of  the  area  surrounding 
event  i,  and  a  number  of  other  proximal  events  not  discussed  in  this  paper, 
clearly  shows  four  open  pit  mines.  All  events  other  than  the  calibration  shots 
were  located  by  Cliff  Thurber  (Thurber  et  al,  1988)  using  the  method  described 
by  Bratt  and  Bache  (1988). 

All  spectral  estimates  displayed  in  this  paper  have  been  computed  using  an 

adaptive  multi-taper  algorithm  described  by  Thompson  (1982).  The  spectral 
leakage  that  plagues  all  estimates,  especially  the  highly-colored  spectra  which 

we  encountered  in  this  study,  is  dealt  with  in  an  effective  manner  by 

employing  several  tapers  that  minimize  the  leakage  from  outside  a  pre¬ 
specified  bandwidth.  Secondly,  the  analysis  of  the  time  dependence  of  the 

energy  in  the  coda  requires  the  computation  of  spectra  from  short  time  series 
which,  especially  at  the  times  of  onsets,  are  particularly  non-stationary. 
Park,  Lindberg  and  Vernon  (1987)  demonstrated  the  multi-taper  method  is 
superior  in  this  situation  to  conventional  single  taper  algorithms  because  it 

employs  several  tapers  that  do  not  overemphasize  the  data  at  the  center  of  the 
time  series  and  discard  significantly  less  data  than  conventional  single-taper 

algorithms. 


DATA  ANALYSIS 

Since  we  are  looking  for  the  presence  or  absence  of  time-independent 
banding,  we  have  found  it  useful  to  calculate  frequency-time  displays  of 
seismograms,  known  as  sonograms  in  acoustics  (Markel  and  Gray,  1976), 
because  they  preserve  the  time  dependence  information. 

It  is  logical  to  begin  the  analysis  with  an  event  about  which  we  know  a  great 
deal.  The  calibration  event  Chemex  2  was  a  20  T  explosion  detonated  at  50.00°  N, 
77.34°  E  at  a  depth  of  17  m  in  a  mining  tunnel  drilled  through  granite  (see 
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table  1).  The  vertical  component  seismogram  for  this  event  recorded  at  high 
gain  by  the  GS-13  seismometer  and  the  corresponding  sonogram  are  displayed 
in  figure  4.  The  velocity  seismogram  has  been  decimated  to  one  point  in  16 
(retaining  the  maximum  and  minimum  amplitudes),  and  is  otherwise 
unfiltered.  The  spectral  estimates  in  this  and  all  subsequent  figures  are  of  the 
acceleration  amplitude  plotted  on  a  linear  scale.  Each  spectral  estimate 
consists  of  a  weighted  sum  of  7  sub-estimates  computed  using  a  time-bandwidth 
product  of  4.  Each  estimate  has  been  computed  from  two  seconds  of  data  (i.e. 
S00  samples),  padded  to  4  seconds  with  zeros.  A  fairly  broadband 

compressional  onset  occurs  27  seconds  after  the  event.  At  this  range  (157  km) 
we  expect  that  the  onset  energy  has  traveled  solely  within  the  crust.  A 
relatively  narrow  band,  high  amplitude  shear  onset  arrives  at  48  seconds 

followed  roughly  10  seconds  later  by  a  small  amplitude  surface  wave  packet. 
The  most  striking  feature  of  this  sonogram  is  that  the  energy  is  clearly  not 
organized  into  time-independent  frequency  bands  but,  with  the  exception  of 
the  onsets,  is  distributed  fairly  randomly. 

We  strongly  suspect  that  event  c  (figure  3)  is  a  ripple-fired  mine  explosion.  It 
was  recorded  by  the  same  seismometer  and  processed  in  exactly  the  same 
manner  as  Chemex  2  and  is  displayed  in  figure  5.  The  sonograms  for  this  event 
and  Chemex  2  are  strikingly  different.  The  energy  is  clearly  organized  into 
time-independent  bands  regularly  spaced  at  roughly  6  Hz.  Considering  the 
regularity  of  the  modulation,  it  seems  reasonable  to  accept,  for  now,  the  initial, 
simplest,  quarry  blast  model  discussed  above.  A  simple  calculation  thus 
suggests  that  this  event  consisted  of  a  number  of  sub-events  spaced  on  the 
order  of  167  ms  from  each  other.  The  actual  number  of  sub-events  involved  is 
not  as  simply  assessed  since  the  estimation  involves  using  the  falloff  rates  of 
the  spectral  peaks  which  are  much  more  difficult  to  compute. 

These  time-independent  bands  are  observed  in  virtually  all  the  unidentified 

events.  The  primary  difference  between  them  is  the  period  of  this  modulation 
along  the  frequency  axis.  A  third  event's  sonogram,  (event  i  in  figure  3, 
located  105  km  from  Bayanaul)  is  displayed  in  figure  6.  The  broader  spectral 
spacing  (df  =  16  Hz)  suggests  a  more  closely  grouped  set  of  sub-shots  separated 
by  roughly  63  ms. 
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In  order  to  enhance  the  time-independent  modulation  we  remove  the  large 
scale  structure  caused  by  the  compressional  and  shear  onsets  as  well  as  the 
high  frequency  spectral  falloff.  We  have  found  that  an  effective  method 

involves  differencing  two  versions  of  each  original  spectral  estimate.  As 
figure  7  illustrates  we  compare  a  relatively  unsmoothed  version  with  one  that 
resolves  only  the  large  scale  structure  in  order  to  extract  the  regular 
modulation.  In  practice,  when  analyzing  the  events  considered  in  this  paper, 
the  smoothed  versions  were  obtained  by  simply  convolving  the  original 

spectral  estimates  with  boxcar  functions  spanning  2.5  and  1.0  Hz.  We 

represent  all  regions  of  the  sonogram  matrix  where  the  local  power  is  high 

relative  to  the  more  regional  average  power  by  a  value  of  +1  (denoted  as  white 
in  this  and  all  subsequent  figures)  and  where  it  is  low  by  a  -1  (black).  In  this 
manner  the  bulk  of  the  magnitude  information  is  discarded  and  the  true 
sonogram  matrix  is  "flattened"  to  a  very  simple,  yet  informative,  binary 
matrix. 

The  binary  matrix  representing  the  first  35  seconds  following  the 

compressional  onset  of  event  Chemex  2  recorded  at  Bayanaul  to  a  frequency  of 
35  Hz  is  displayed  in  Figure  8.  Aside  from  the  suggestion  of  time-independent 
structure  in  the  vicinity  of  5  Hz,  the  energy  is  distributed  more  or  less 
randomly.  It  seems  intuitively  reasonable  that  a  single  source  explosion 
should  give  rise  to  a  coda  in  which  the  energy  is  distributed  randomly  as  a 
function  of  frequency  and  time  if  the  dominant  coda  generation  mechanism 
involves  the  scattering  of  energy  by  a  random  distribution  of  randomly  sized 

inhomogeneities.  This  example  suggests  that  layer  resonance,  proposed  by 
Sereno  and  Orcutt  (1985a)  as  the  dominant  source  of  coda  in  the  oceanic 
lithosphere,  plays  only  a  supporting  role  in  generating  seismic  coda  in  this 
region  of  the  continental  crust.  We  contrast  this  result  with  the  coda  produced 
by  event  c  (figure  9)  which  displays  obvious  time-independent  spectral 

modulation. 

We  examined  a  fourth  event  (event  d  in  Figure  3)  to  find  that  the  regular 
modulation  is  not  restricted  to  the  vertical  component.  Indeed,  as  Figure  10 
reveals,  the  modulation  appears  to  be  quite  independent  of  the  component  of 
displacement.  Only  the  vertical  and  radial  components  are  shown,  but  the 

transverse  component  is  nearly  identical.  This  observation  lends  support  to 
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the  hypothesis  that  the  time-independent  spectral  modulation  is  a  source  and 
not  a  propagation  effect.  This  observation  is  not  unique  to  this  event  but 
appears  to  be  shared  by  all  events  other  than  the  calibration  shots. 

The  quarry  blast  spectral  modulation  is  not  only  independent  of  component 
but  appears  to  be  independent  of  the  azimuth  from  the  source  to  the  receiver 
as  well  as  source- receiver  range.  In  figure  11  we  compare  the  vertical 
component  modulation  of  event  d  as  recorded  at  Bayanaul  (at  a  range  of  159  km 
and  a  back  azimuth  of  268°)  and  Karkaralinsk  (at  a  range  of  217  km  and  a  back 
azimuth  of  317°).  The  independence  of  azimuth  suggests  that  the  array  of  sub¬ 
shots  comprising  event  d  may  be  very  small.  The  apparent  time  offsets  caused 
by  spacing  the  sub-shots  become  significant  when  the  array  is  large  and  as 
discussed  earlier  these  apparent  offsets  have  an  azimuthal  dependence. 

To  constrain  the  probable  dimensions  of  a  typical  quarry  blast  we  employ  the 
formula  of  Smith  (1989)  to  equate  true  spatial  offsets  with  apparent  time 
offsets  8Tja 

ST3  =  p  ,  /  5X2sin20  +  5Y2cos20  +  5T.  (7) 

i  r  v  i  t  t 


Each  sub-event  occurs  at  point  (8Xj,8Yj)  in  space  and  at  time  8Tj.  The  energy 
under  consideration  travels  at  a  slowness  of  p  s  km'*  along  an  azimuth  from 
the  source  to  the  receiver  of  0  degrees.  To  permit  a  crude  calculation  we  adopt 
an  array  consisting  of  just  two  shots.  We  simplify  the  calculation  by  aligning 
the  X  and  Y  axes  parallel  and  perpendicular  to  the  line  joining  the  shots 
respectively.  To  use  this  equation,  we  must  first  estimate  the  dependence  of 

the  modulation  pattern  on  source  to  receiver  azimuth.  In  general  it  is  true  in 
this  dataset  that  little  difference  in  the  modulation  period  associated  with  any 

single  event  can  be  discerned  between  the  recordings  made  at  Karkaralinsk 
and  Bayanaul.  Thus,  for  a  crude  estimate  of  the  likely  source  dimensions  we 
model  a  typical  event  (d,  figure  10)  with  sub-shots  arranged  spatially  to  give 
rise  to  the  maximum  possible  azimuthal  dependence  of  the  spectral  modulation 

when  observed  at  these  two  stations.  We  will  realize  the  maximum  possible 
change  in  modulation  period  between  Bayanaul  and  Karkaralinsk  if  we  place 
this  dipolar  shot  array  so  that  the  normal  to  the  line  joining  the  shots  lies 
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exactly  half  way  between  the  two  receivers.  Using  elementary  error  analysis, 
and  knowing  that  8T=1/8F  (where  8T  is  the  time  offset  and  8F  is  the  spacing  of 

adjacent  peaks  in  the  frequency  domain)  we  Find  that  the  change  in  8T  (or 

2  2 
8  T)  required  by  a  change  in  the  modulation  spacing  (8  F)  is  given  by: 

S2T  =  -  (8) 

(sf) 

If  we  observe  n  cycles  in  the  frequency  span  of  AF  Hz  in  one  recording  of  an 
event,  and  in  AF  ±  dF  Hz  in  another  recording  then 

82T  =  -  (9) 

(AF) 

Combining  equations  (7)  and  (9)  we  find  that 

SX  - —  r — adF - T  (10) 

(AF)  |  sin(0  )  -  sin(6  )  j 
v  L  1  z  J 

where  8X  is  the  inferred  shot  spacing  and  01  and  02  are  the  azimuths  to 

Bayanaul  and  Karkaralinsk  (and  equal  +24.6°  and  -24.6°  respectively).  In 
figure  11  we  observe  11  cycles  present  in  a  frequency  band  of  30  ±  1  Hz  in  both 
the  recordings  of  this  event.  Thus,  using  equation  (10),  and  considering  a 

typical  crustal  ray  with  a  slowness  of  1/7  s  km'*  we  infer  a  shot  spacing  of  100 
m  or  less.  Clearly  many  assumptions  have  been  made  in  this  calculation  so  it 
should  be  regarded  as  a  crude  estimate  at  best.  We  have  no  a-priori 

information  on  what  blasting  techniques  (shot  geometry  and  timing)  the 

Soviet  mining  engineers  employ  in  their  quarry  blasts.  Combining  this  fact 

with  the  observations  that  the  spectral  modulations  are  regular  in  frequency 
and  independent  of  time  and  source-receiver  azimuth  it  seems  inappropriate 
to  attempt  a  more  sophisticated  modeling  by  allowing  spatial  offsets.  We  can 
adequately  model  these  Soviet  quarry  blasts  by  employing  either  of  the  two 
models  proposed  earlier  in  this  paper. 

Assuming  linear  superposition  (and  the  first  quarry  blast  model),  it  is  a  simple 
matter  to  synthesize  a  quarry  blast  by  linearly  superposing  a  Green's  function 
upon  itself  after  offsetting  by  a  delay  time  prior  to  each  addition.  By  chance  it 
turns  out  that  Chemex  2  is  at  nearly  the  same  distance  from  Bayanaul  as  is 
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event  d,  and  thus  it  is  eligible  for  use  as  a  "Green's  function"  if  lateral 
variations  of  crustal  structure  are  small.  This  observed  rather  than 
synthesized  Green's  function  has  the  obvious  drawback  of  being  produced  by  a 
significantly  larger  source  than  a  likely  sub-event  in  a  typical  quarry  blast, 
so  that  its  corner  frequency  will  be  lower.  This  shouldn't  be  a  significant 
problem  since,  as  Evernden,  Archambeau  and  Cranswick  (1986)  note,  the 
Sharpe  explosion  model  (Sharpe,  1942),  predicts  the  corner  frequency  of  a 
tamped  explosion  of  this  size  to  be  30  to  35  Hz.  This  is  approximately  the  upper 
frequency  limit  in  ail  binary  sonograms  considered  in  this  study.  The 
usefulness  in  general  of  such  a  Green's  function  is  limited  because  it  is  only 
available  at  a  small  number  of  ranges  (defined  by  the  offsets  between  the  two 
events  and  the  two  receivers).  It  has  the  distinct  advantage  vis-a-vis 
synthetics  that  the  propagation  transfer  function  is  nearly  exact  (assuming 
horizontal  layers).  By  taking  the  vertical  component  of  motion  from  Chemex  2 
observed  at  Bayanaul  (figure  4)  and  stacking  it  upon  itself  twice  after 

offsetting  each  successive  trace  by  133  ms,  we  obtain  the  "synthetic"  quarry 
blast  displayed  in  figure  12.  The  binary  sonogram  corresponding  to  a 
synthetic  quarry  blast  consisting  of  3  shots  offset  by  382  ms,  displayed  in 
figure  13,  compares  well  with  the  actual  binary  sonogram  computed  from 
event  d.  Although  non-linear  effects  must  clearly  occur  during  a  quarry 
blast,  simple  linear  theory  can  be  very  effective  in  reproducing  the  observed 

modulation.  In  principle,  we  can  time  an  observed  ripple  fired  quarry  blast 

repetition  rate  with  this  method.  In  this  case  we  have  illustrated  a  likely  offset 
time  (382  ms)  for  event  d. 

DISCRIMINATION  BETWEEN  SINGLE  AND  MULTIPLE  SHOT  EXPLOSIONS 

The  method  described  above  of  reducing  sonograms  to  binary  matrices 
provides  simple  patterns  that  allow  visual  discrimination  between  quarry 
blasts  and  simpler  events  (single  event  explosions  and  perhaps  earthquakes). 
Clearly  there  are  a  number  of  mathematical  operations  we  can  employ  to 
reduce  these  matrices  to  scalars  which  reflect  the  presence  or  absence  of 

time-independent  spectral  banding.  Two  methods  seem  promising  at  this  time 
but  have  not  yet  been  developed.  One  involves  the  comparison  of  three 
component  binary  matrices  by  performing  a  3  way  cross-correlation.  As 
discussed  above  the  quarry  blast  modulation  is  a  source  effect  and  is  not 
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dependent  on  the  recorded  component.  We  expect  a  relatively  high  cross¬ 
correlation  between  the  three  components  of  recorded  motion  resulting  from 
a  quarry  blast  relative  to  that  computed  from  simpler  events.  In  an  alternate 
method  the  2-dimensional  Fourier  transform  of  the  binary  sonogram  matrix  is 
computed.  We  expect  the  presence  or  absence  of  time-independent  modulation 
would  be  reflected  in  the  time  axis  DC  power  levels,  those  events  with  time- 
independent  modulation  should  have  relatively  high  levels.  For  the  present 

purposes,  however,  we  feel  it  is  sufficient  to  recognize  simply  that  distinct 
time  versus  frequency  pattern  differences  between  quarry  blasts  and 
calibration  explosions  exist.  By  comparing  individual  columns  in  the  binary 
sonogram  matrices  calculated  from  single  and  multiple-event  explosions  we 
can  assess  the  enhanced  discrimination  potential  of  binary  sonogram  matrices 

over  that  of  single  spectral  estimates.  In  figure  8  many  of  the  columns  possess 
a  regular  modulation  comparable  to  that  seen  in  single  columns  computed 
from  the  coda  resulting  from  a  quarry  blast  (for  example  figure  9).  Individual 
spectral  estimates  computed  from  even  simple  single  event  explosions  can 
have  significant  regular  modulation  and  thus  distinguishing  between  these 
and  more  complex  multiple-event  explosions  using  single  spectra  is  expected 
to  be  difficult.  The  profound  differences  between  these  two  types  of  events 
becomes  obvious  when  many  spectra  offset  from  each  other  in  time  are 
computed.  The  persistence  of  spectral  features  through  time  in  the  coda  of 

multiple  source  events,  but  not  single  source  events,  is  what  sets  them  apart. 

A  weakness  of  this  discriminant  and  indeed  any  algorithm  which  seeks  to 
separate  quarry  blasts  from  other  events  on  the  basis  of  their  unusual  spectral 
color  is  that  time-independent  spectral  modulation  can  be  acquired  during 
propagation.  The  binary  sonogram  matrix  of  Chemex  2  recorded  at  Bayanaul 

possesses  a  faint  but  undeniable  time-independent  structure  at  about  5  Hertz. 
Comparison  of  this  figure  with  the  binary  sonogram  computed  from  Chemex  2 

recorded  at  Karkaralinsk  reveals  that  the  two  stations  have  recorded 

essentially  the  same  patterns.  The  same  exercise  performed  on  the  Chemex  1 
recordings  at  Bayanaul  and  Karkaralinsk  revealed  faint  but  strikingly  similar 
time-independent  patterns  distinct  from  the  Chemex  2  patterns.  These 

features  are  most  likely  a  source  effect  since  similar  patterns  have  sources, 

not  receivers,  in  common.  Since  the  Chemex  events  are  known  to  be  single 
source  events,  the  most  plausible  explanation  is  that  the  patterns  are  due  to 
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reverberations  occurring  in  a  shallow,  near  source  horizon.  Although  these 
time-independent  effects  are  faint  they  are  worth  considering  further. 
Reverberations  can  occur  at  any  time  during  propagation  from  the  source  to 

the  receiver,  but  the  two  points  at  which  the  effects  are  going  to  be  most 

noticeable  and  most  likely  is  near  the  source  and  near  the  receiver.  These  are 
the  only  two  points  at  which  the  recorded  body  wave  energy  is  most  likely  to 
encounter  a  low  velocity  horizon.  Prior  to  using  the  discriminant,  we  are 

proposing  that  a  receiver  must  be  shown  to  be  free  of  reverberations  by 
recording  single  event  explosions.  If  the  site  is  able  to  record  these  events 

with  no  time-independent  effects  then  it  is  eligible  for  this  type  of  work. 
Fortunately  the  Chemex  events  have  not  produced  significant  time- 
independent  structure  (the  high  frequency  overtones  have  not  developed), 

but  they  do  indicate  that  it  may  be  possible  to  disguise  a  single  event  explosion 
as  a  quarry  blast  by  detonating  it  in  a  low  velocity  horizon. 

A  second  problem  exists  when  events  have  been  recorded  at  a  time  of  high 
noise  levels,  especially  when  the  events  are  small  and/or  distant  and  when 
spectral  spikes  exist  in  the  noise.  Prominent  features  in  ground  or 
instrumental  noise  can  emerge  from  the  signal  after  the  onset  and  be 
misinterpreted  as  time-independent  features  due  to  the  source,  especially  if 

more  than  one  harmonic  is  present.  To  combat  this  problem  it  is  necessary  to 
obtain  and  analyze  pre-event  noise  samples  to  permit  the  identification  of 
time-independent  features  in  the  coda  that  are  really  due  to  noise.  In  this 
study,  all  time-independent  features  attributed  to  the  source  are  seen  on  both 
the  Karkaralinsk  and  Bayanaul  records.  It  is  also  true  that  while  signal 

spectra  between  the  two  receivers  have  high  coherence,  the  pre-event  noise 
spectra  do  not. 

In  practice  it  is  clearly  desirable  that  more  than  one  station  be  used  when 
attempting  to  discriminate  quarry  blasts  from  simpler  events.  As  discussed 
earlier,  if  the  quarry  blast  sub-shot  array  is  not  infinitely  small  we  expect  the 
modulation  pattern  to  show  an  azimuthal  dependence.  There  are  certain  array 
geometries,  for  example  linear  arrays,  that  will  produce,  at  certain  source- 
receiver  azimuths,  negligible  apparent  time  offsets  between  the  sub-shots  and 
thus  little  or  no  time-independent  spectral  structure.  It  is  extremely  unlikely 

that  a  second  receiver  simultaneously  recording  this  event  at  a  different 
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azimuth  will  also  have  such  a  problem.  None  of  the  quarry  blasts  we  have 
examined  show  evidence  for  such  azimuthal  effects. 


CONCLUSIONS 

We  have  observed  a  significant  difference  between  events  which  consist  of  a 
number  of  sub-explosions  closely  grouped  in  space  and  time  and  simpler, 
single  explosion,  events.  The  former  type  of  event  (a  quarry  blast)  possesses  a 
distinctive  time-independent  spectral  structure  that  can  be  easily  explained 
assuming  linear  superposition  of  the  wavefields  produced  by  each  of  the  sub¬ 
events.  The  latter  type  of  source  does  not  cause  a  regular  organization  of  the 
energy  and  thus  the  distribution  of  energy  in  the  coda  is  controlled  by  the 
propagating  medium.  The  observed  randomness  of  this  energy  as  a  function 
of  frequency  and  time  suggests  that  the  dominant  mechanism  of  coda 
generation  in  this  region  of  the  continental  crust  involves  scattering  by  small 
scale  inhomogeneities.  We  conclude  these  two  types  of  events  can  be  most 
effectively  distinguished  from  each  other  by  exploiting  the  expected 
differences  in  the  time-  and  frequency-evolution  of  energy  in  the  onsets  and 
coda.  We  propose  a  time  versus  frequency  pattern-based  discriminant. 

Although  we  have  no  a-priori  information  on  the  mining  practice  employed 
by  the  Soviet  quarry  blast  engineers  we  infer  regular  offsets  in  time  and  small 
spatial  array  dimensions  by  observing  that  the  modulations  are  regular  in 
frequency  and  independent  of  time  and  source-receiver  azimuth. 

In  the  future  we  will  expand  this  analysis  to  include  earthquakes  and 
explosions  in  other  regions.  We  expect  we  can  demonstrate  the  independence 
of  the  algorithm  on  mining  practice  and  geologic  setting.  We  hope  to  extend 
the  pattern-based  discriminant  by  reducing  the  binary  sonogram  matrices  to 
scalars  which  reflect  the  presence  or  absence  of  time-independent 
modulation. 
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FIGURE  CAPTIONS 


Figure  1.  Spectral  modulations  predicted  for  events  consisting  of 

1  to  5  sub-shots  spaced  evenly  at  100  ms. 

Figure  2.  Spectral  modulation  predicted  for  an  event  consisting 

of  SO  sub-shots.  The  offset  times  are  distributed  as  a 
Gaussian  with  a  mean  of  63  ms  and  a  variance  of  6.3 


Figure  3.  Receiver  and  event  Geography. 

Figure  4.  Seismogram  resulting  from  chemex  2  recorded  at  a 

range  of  157  km  by  the  vertical  surface  seismometer 
at  Bayanaul  and  corresponding  sonogram. 

Figure  5.  Seismogram  resulting  from  event  c  recorded  at  a  range 

of  264  km  by  the  vertical  surface  seismometer  at 
Bayanaul  and  corresponding  sonogram. 

Figure  6.  Seismogram  resulting  from  event  i  recorded  at  a  range 

of  105  km  by  the  vertical  surface  seismometer  at 
Bayanaul  and  corresponding  sonogram. 

Figure  7.  Original  spectral  estimate  and  two  versions,  one 

relatively  unsmoothed  and  the  other  relatively  highly 
smoothed.  This  figure  is  intended  to  illustrate  the 
means  by  which  we  reduce  each  spectral  estimate  to  a 
"binary  spectral  estimate".  The  regions  of  locally  high 
power  are  represented  by  a  +1  (white).  Regions 
deficient  in  power  are  represented  by  a  -1  (black). 

Figure  8.  Vertical  component  binary  sonogram  matrix 

corresponding  to  chemex  2  recorded  at  Bayanaul.  The 
first  25  seccrds  of  coda  after  the  compressional  onset 
are  represented  in  this  and  all  subsequent  figures. 

Figure  9.  Vertical  component  binary  sonogram  matrix 

corresponding  to  event  c  recorded  at  Bayanaul. 

Figure  10.  Binary  sonogram  matrices  representing  the  vertical 

(a)  and  radial  (b)  components  of  displacement 
recorded  at  Bayanaul  from  event  d. 

Figure  11.  Vertical  component  binary  sonogram  matrices 

corresponding  to  event  d  recorded  at  Bayanaul  (a)  and 
Karkaralinsk  (b). 

Figure  12.  (a)  Time  series  and  sonogram  corresponding  to 

chemex  2  recorded  at  Bayanaul.  (b)  Time  series  and 
sonogram  representing  a  "synthetic"  quarry  blast 
involving  3  shots  spaced  in  time  by  133  ms. 
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Figure  13.  (a)  Vertical  component  binary  sonogram  matrix 

corresponding  to  a  synthetic  quarry  blast  consisting  of 
3  sub-events  spaced  in  time  by  382  ms.  (b)  Vertical 
component  binary  sonogram  matrix  corresponding  to 
event  d  recorded  at  Bayanaul. 
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Table  1.  Event  Locations 


Event 

Latitude 

Longitude 

Origin  Time 

(°N) 

(°B) 

year/day/h:m:s 

CH2 

50.000 

77.340 

87/245/09:27:04.950 

EVC 

49.333 

72.667 

87/135/10:35:00.2 

EVD 

50.742 

73.275 

87/141/09:16:43.3 

EVI 

51.760 

75.572 

87/146/08:33:26.5 
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Figure  10 
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Figure  11 
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Figure  13 
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OBJECTIVE 


The  primary  objective  of  this  study  is  to  develop  an  algorithm  that  will 
routinely  discriminate  between  ripple-fired  mine  explosions  (events  that 
consist  of  several  sub-explosions  closely  grouped  in  space  and  time)  and  all 
other  seismic  events  (including  earthquakes  and  single-event  explosions).  To 
develop  the  algorithm  and  test  its  independence  from  local  mining  practice 
and  geologic  setting  we  have  selected  two  independent  data  sets  (data  collected 
by  the  NORESS  array  in  southern  Norway  and  by  sensors  deployed  in 
Kazakhstan  U.S.S.R.).  To  assess  the  relative  value  of  these  two  data  sets  for  this 
type  of  work  we  have  conducted  a  comparative  analysis  of  ambient  noise  levels 
at  these  two  sites. 
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SUMMARY 


A  commonly  used  technique  in  mining  practice,  known  as  ripple-firing, 
involves  the  detonation  of  a  number  of  sub-explosions  offset  from  each  other 
by  small  distances  and  times.  The  intent  of  this  procedure  is  to  enhance 
fracturing  of  the  rock  and  reduce  ground  motion  in  areas  proximal  to  the 
mine.  The  observation  has  been  made  by  several  researchers  (Smith,  1989, 
Stump  and  Reamer,  1988,  Baumgardt  and  Ziegler,  1988)  that  such  events  tend  to 
give  rise  to  seismic  coda  possessing  highly  colored  spectra.  That  is,  spectra 
enriched  in  power  at  certain  preferred  frequencies  and  depleted  in  power  at 

other  frequencies. 

This  spectral  color  is  certainly  due  to  the  interaction  of  the  time-offset 
wavefields  produced  by  each  sub-explosion.  The  manner  in  which  the 
wavefields  interact  is  not  known,  and  undoubtedly  is  largely  non-linear, 
however  we  feel  that  simple  linear  theory  is  sufficient  to  describe  the  most 

obvious  result,  specifically  the  pronounced  spectral  modulation.  Briefly,  the 
regular  repetition  and  superposition  of  similar  seismic  motions  in  the  time 

domain  leads  to  regular  amplification  and  suppression  of  power  in  the 
frequency  domain.  Furthermore,  as  is  demonstrated  by  Hedlin  et  al  (1989)  in 

the  limit  where  all  sub-explosions  occur  at  the  same  point  in  space,  the 

spectral  features  are  independent  of  time  in  the  coda.  Because  we  predict  that 
ripple-fired  events  should  give  rise  to  pronounced,  time-independent,  spectral 
features,  we  have  found  it  useful  to  calculate  frequency-time  displays  of  the 

seismograms,  known  as  sonograms. 

In  figure  1  we  display  the  vertical  component  time  series,  and  sonogram,  of  a 

20  T  single-event  calibration  explosion  detonated  and  recorded  in  Kazakhstan 
U.S.S.R.  in  1987.  The  recording  was  made  by  a  GS-13  seismometer  sampled  at 
250  s*1  and  deployed  at  the  surface  near  the  town  of  Bayanaul.  The  energy  in 
the  coda  from  this  event  is  distributed  fairly  randomly  as  a  function  of 

frequency  and  time.  There  is  no  evident  organization  of  the  energy  into  time- 
independent  bands.  In  figure  2  we  show  the  result  of  the  same  processing 
techniques  applied  to  an  unidentified  event  that  occurred  264  km  from  the 
station.  The  spectra  of  this  event  are  distinctly  modulated  and  furthermore  the 
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energy  is  very  clearly  organized  into  time-independent  bands.  We  strongly 
suspect  that  this  event  was  a  ripple-fired  quarry  blast. 

We  have  found  it  useful  to  suppress  the  large  scale  structure  provided  by  the 
compressional  and  shear  onsets  as  well  as  the  high  frequency  spectral  falloff 
and  thus  enhance  any  time-independent  features  that  may  be  present  by 
comparing  two  versions  of  each  spectral  estimate.  As  figure  3  illustrates  we 
compare  a  relatively  unsmoothed  version  of  the  original  spectral  estimate 
with  one  that  is  heavily  smoothed.  We  represent  all  regions  of  the  sonogram 
matrix  where  the  local  power  is  high  relative  to  the  regional  power  by  a  value 
of  +1  (denoted  as  white  for  display  purposes)  and  where  it  is  low  by  a  -1 
(black).  In  this  manner  we  discard  the  bulk  of  the  magnitude  information  and 
"flatten"  the  original  sonogram  matrix  to  a  very  simple,  yet  informative, 
binary  matrix.  In  figures  4  and  S  we  display  the  binary  sonogram  matrices 

corresponding  to  the  originals  displayed  in  figures  1  and  2  respectively.  This 
method  of  reducing  sonogram  matrices  to  binary  form  provides  simple 
patterns  that  allow  visual  discrimination  between  mine  explosions  and  single¬ 
event  explosions. 

To  extend  this  work  we  have  examined  data  collected  by  the  NORESS  array  in 
southern  Norway.  Our  goal  was  to  discover  whether  earthquakes  can  be  as 
easily  distinguished  from  mine  blasts  as  single-event  explosions.  We  also 
wanted  to  determine  if  the  algorithm  would  be  successful  in  a  different 
geologic  setting  and  in  an  area  of  different  mining  practice  than  present  in 

central  Asia.  We  have  concentrated  our  efforts  thus  far  on  events  known  to 

either  be  earthquakes  or  mine  explosions.  The  NORESS  array  data  have  a 
distinct  advantage  over  the  Kazakh  data  since  typically  30  sensors 
simultaneously  record  each  event.  The  obvious  disadvantage  is  that  the  array 
data  are  only  digitized  at  40  s‘*  and  thus  cannot  resolve  spectral  features  above 
20  Hz.  In  figure  6  we  display  the  sonogram  of  a  presumed  earthquake  which 
occurred  454  km  from  the  array.  This  display  has  been  produced  by  stacking 
the  amplitude  spectra  from  each  of  the  24  vertical  channels  after  converting 

each  spectrum  to  binary  form  and  assuming  infinite  phase  velocity.  The 
compressional  onset  occurs  at  12  seconds.  As  expected,  little  structure  can  be 
discerned  other  than  the  microseismic  noise  below  2  Hz.  Figure  7  was 
computed  from  the  coda  produced  by  a  Blasjo  mine  explosion  which  occurred 
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302  km  from  the  array.  The  compressional  onset  occurs  at  15  seconds  and  is 
followed  by  a  very  obvious  time-independent  banding.  The  time-independent 
structure  present  prior  to  the  onset  is  due  to  long  lived  spectral  lines  in  the 
noise  and  does  not  correlate  with  the  spectral  features  present  in  the  coda. 

In  figures  8  and  9  we  display  the  ambient  noise  levels  present  at  the  high- 
frequency  element  of  the  NORESS  array  (site  NRFO  ;  sampling  at  125  s'*)  and  in 
the  boreholes  at  the  Kazakh  sites  near  the  towns  of  Bayanaul  and 
Karkaralinsk.  All  seismometers  recorded  the  vertical  component  of  velocity 
and  were  deployed  at  60,  66  and  99  meters  depth  respectively.  The  noise  levels 
are  represented  by  three  curves  which  illustrate  the  average  and  95% 
confidence  limits  of  acceleration  power.  Above  1  to  2  Hz  the  Kazakh  noise 
levels  are  generally  slightly  (up  to  10  dB)  higher.  The  apparent  rise  in  power 
above  40  Hz  at  the  NORESS  array  is  probably  not  natural  but  due  to  digitization 
noise.  Below  1  to  2  Hz  the  NORESS  noise  levels  are  significantly  higher  than  in 
central  Asia.  Although  these  curves  do  not  extend  across  the  entire 
microseism  band,  the  suggestion  is  that  microseismic  noise  levels  are  higher 
in  Norway  than  in  central  Asia.  This  seems  intuitively  reasonable  given  the 
relative  proximity  of  the  NORESS  array  to  an  ocean.  This  appears  to  be  a  long- 
range  extension  of  the  observation  made  by  Hedlin  et  al  (1988)  that 
microseismic  noise  levels  appeared  to  be  inversely  correlated  with  the 
distance  from  the  shoreline  in  the  range  of  1  to  100  km. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


We  observe  significant  differences  between  ripple-fired  mine-explosions  and 
earthquakes  and  single-event  explosions.  The  close  grouping  in  space  and 
time  of  several  sub-explosions  appears  to  be  responsible  for  the  impression  of 
prominent  time-independent  spectral  modulation  on  the  frequency  spectrum 
of  the  resultant  coda. 

We  plan  to  extend  further  the  pattern-based  discriminant  by  reducing  the 
binary  sonogram  matrices  to  scalars  which  reflect  the  presence  or  absence  of 
time-independent  modulation. 

It  appears  that  lower  microseismic  noise  levels  will  be  realized  at  sites  far 
removed  from  large  bodies  of  water. 
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FIGURE  CAPTIONS 


Figure  1.  Seismogram  resulting  from  a  single-event  calibration  explosion 
recorded  at  a  range  of  157  km  by  the  vertical  surface  seismometer  at  Bayanaul 
and  corresponding  sonogram. 

Figure  2.  Seismogram  resulting  from  a  suspected  mine-explosion  recorded 
at  a  range  of  264  km  by  the  vertical  surface  seismometer  at  Bayanaul  and 
corresponding  sonogram. 

Figure  3.  Original  spectral  estimate  and  two  versions,  one  relatively 
unsmoothed  and  the  other  relatively  highly  r  oothed.  This  figure  is  intended 
to  illustrate  the  means  by  which  we  reduce  each  spectral  estimate  to  a  "binary 
spectral  estimate".  The  regions  of  locally  high  power  are  represented  by  a  +1 
(white).  Regions  deficient  in  power  are  represented  by  a  -1  (black). 

Figure  4.  Vertical  component  binary  sonogram  matrix  corresponding  to 
the  single-event  calibration  explosion  recorded  at  Bayanaul  (refer  to  figure 
1).  The  first  35  seconds  of  coda  after  the  compressional  onset  are  represented 
in  this  and  the  following  figure. 

Figure  5.  Vertical  component  binary  sonogram  matrix  corresponding  to 
the  suspected  mine-explosion  recorded  at  Bayanaul. 

Figure  6.  Vertical  component  sonogram  matrix  corresponding  to  a 
presumed  earthquake  recorded  at  the  NORESS  array. 

Figure  7,  Vertical  component  sonogram  matrix  corresponding  to  a  Blasjo 
mine  explosion  recorded  at  the  NORESS  array. 

Figure  8.  The  average  and  95%  confidence  limits  of  ambient  noise  levels 
recorded  at  the  NORESS  array  and  near  the  town  of  Bayanaul  in  Kazakhstan, 
U.S.S.R. 

Figure  9.  The  average  and  95%  confidence  limits  of  ambient  noise  levels 
recorded  at  the  NORESS  array  and  near  the  town  of  Karkaralinsk  in 
Kazakhstan,  U.S.S.R. 
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Figure  3.  Original  spectral  estimate  and  two  versions,  one  relatively  unsmoothed  and 
the  other  relatively  highly  smoothed.  This  figure  is  intended  to  illustrate  the  means  by 
which  we  reduce  each  spectral  estimate  to  a  "binary  spectral  estimate".  The  regions  of 
locally  high  power  are  represented  by  a  +1  (white).  Regions  deficient  in  power  are 
represented  by  a  -1  (black). 
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Fi{ur«  4.  Vertical  component  binary  sonogram  matrix  corresponding  to  the  single-  Fijure  5.  Vertical  component  binary  sonogram  matrix  corresponding  to  the  suspected 

event  calibration  explosion  recorded  at  Bayanaul  (refer  to  figure  1).  The  first  35  seconds  mine-explosion  recorded  at  Bayanaul. 

of  coda  after  the  cor-otessional  onset  are  represented  in  this  and  the  following  figure. 
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ABSTRACT 


We  observe  significant  modulation  in  spectra  computed  from  seismic  coda 
recorded  by  seafloor  seismometers  and  hydrophones  deployed  in  the  S.W. 
Pacific  and  in  the  Wake  Island  Hydrophone  array.  These  modulations  are 
consistent  with  those  present  in  synthetic  spectra,  generated  using  an  Earth 
model  consisting  of  homogeneous  layers,  and  consistent  with  the  locations  of 
spectral  peaks  we  predict  using  simple  layer  resonance  theory.  Our 
observations  agree  with  those  made  by  Sereno  and  Orcutt  (1985  a.b)  and  are 
consistent  with  their  model  of  coda  generation  by  layer  resonance. 
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INTRODUCTION 


For  years  researchers  have  recorded  seismic  energy  that  has  propagated 
through  the  oceanic  lithosphere.  Distinguishing  features  of  this  energy  are 

that  high  frequencies  are  retained  at  teleseismic  ranges  from  the  source 
(Walker  et  al,  1978).  In  addition  long  coda  are  observed  following  the 

compressional  and  shear  onsets  (Walker  et  al,  1983).  It  is  widely  accepted  that 
this  energy  remains  within  the  oceanic  lithosphere  (Sereno,  1986);  however, 
the  mechanism  by  which  the  seismic  coda  is  generated  has  not  been 
established.  Richards  and  Menke  (1983)  and  Menke  and  Chen  (1984)  propose 
that  scattering  by  small  scale  inhomogeneities  is  an  effective  mechanism  for 

the  generation  of  long  coda.  Recently  an  alternative  theory  was  proposed  by 
Sereno  and  Orcutt  (1985a, b).  They  suggest  that  resonance  of  energy  in  the 

water  column  and  underlying,  low  velocity,  sediment  layer  is  sufficient  to 
generate  the  long  coda  observed.  A  likely  feature  of  coda  produced  by  such  a 
mechanism  is  that  is  should  possess  a  spectral  scalloping,  or  regular 
modulation  in  the  frequency  domain,  produced  by  the  superposition  of 
resonant  phases  in  time.  The  objective  of  this  paper  is  to  reproduce  the  results 
of  Sereno  and  Orcutt  (1985a,b)  by  observing  this  spectral  scalloping. 


RESONANCE  THEORY 


The  theory  behind  the  generation  of  seismic  coda  by  layer  resonance  was 
developed  by  Sereno  and  Orcutt  (1985b)  and  thus  only  a  terse  description  is 
given  here.  The  authors  demonstrate  the  significant  influence  of  low  velocity 
horizons  in  the  oceanic  environment  on  coda  duration  and  amplitude.  This 
layer  resonance  enriches  the  spectrum  of  the  coda  at  certain  equispaced 
eigenfrequencies,  determined  by  the  layer  thickness  and  velocity. 
Considering  a  relatively  simple  case  in  which  a  wavelet  w(t)  is  reverberating 
in  a  single  low  velocity  surface  layer,  the  linear  sum  of  all  reverberations  x(t) 
equals; 


x(t)  =  w(t)  * 


•  e-“"' 


H(t)  •  e 


-int/T 


(1) 


The  time  between  successive  reverberations  is  represented  by  T  and,  to  a  good 
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approximation,  equals  the  two-way  vertical  travel  time  in  the  horizon.  The 
"shah”  fv  iction  is  tapered  exponentially  to  simulate  the  effect  of  attenuation 
and  is  multiplied  by  the  Heaviside  step  function  to  exclude  any  acausal  energy. 
The  final  exponential  is  required  by  the  n  phase  shift  incurred  by  reflection  at 
the  free  surface. 


By  Fourier  transforming  the  preceding  equation  the  authors  found  that  the 
spectrum  of  the  original  wavelet  is  multiplied  by  an  infinite  set  of  staggered 
tapering  functions  which  decay  at  a  rate  controlled  by  and  proportional  to  the 
attenuation  parameter  a.  The  amplitude  spectrum  of  the  sum  of  wavelets  is 
given  by: 


IX(f)l  =  IW(f)l 


III(fT  +  * 


a  -  i2nf 
(a2  +  47t2  f2) 


(2) 


The  thicker  -  or  equivalently,  the  slower  -  the  horizon  the  more  closely  spaced 
the  spectral  highs  will  be. 


In  the  Oceanic  setting  there  are  two  low  velocity  strata,  the  water  column  and 
the  underlying  sediment  layer.  It  is  resonance  in  these  horizons  that  Sereno 
and  Orcutt  (1985a. b)  invoked  in  their  explanation  of  the  origin  of  the 
significant  Pn  and  Sn  coda.  It  is  modulations  from  these  resonances  that  we 
seek  in  the  real  and  synthetic  data  sets  considered  below  as  we  attempt  to 
confirm  their  conclusions. 

DATA  SETS 

The  model  proposed  by  Sereno  and  Orcutt  (1985a, b)  to  explain  the  propagation 

of  the  Pn  phase,  and  the  development  of  long  coda,  was  established  primarily 
using  data  collected  by  Ocean  Bottom  Seismometers  (OBS’s)  and  Hydrophones. 

This  data  set  was  collected  during  Scripps’  Ngendei  expedition  in  the  South- 
West  Pacific  in  1983  (Orcutt  et  al,  1987).  In  this  paper  we  compare  the  results 

of  this  experiment  with  those  of  an  independent  data  set  collected  by  near¬ 

seafloor  sensors  in  the  Wake  Island  Hydrophone  (W1H)  array  (McCreery, 
1987).  We  have  obtained  recordings  of  70  teleseismic  events  made  by  this 
array.  One  advantage  of  this  data  set  is  that  each  event  is  simultaneously 
recorded  by  several  (usually  5)  seafloor  sensors  and  thus  we  can  gain  a  limited 
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improvement  in  signal  to  noise  ratios  by  stacking  in  the  time-domain.  The 
slight  disadvantage  is  that  the  sensors  are  deployed  in  the  water  column.  We 
investigate  the  feasibility  of  testing  a  crustal  propagation  model  using  sensors 
uncoupled  from  the  seafloor  and  we  compare  synthetic  seismograms  computed 
with  receivers  in  both  the  water  column  and  the  underlying  sediment  layer. 
We  have  used  the  Wavenumber  Integration  algorithm  for  this  purpose  (Luco 
and  Apsel,  1983;  Apsel  and  Luco,  1983).  This  algorithm  requires  that  the  Earth 
be  modeled  as  a  stack  of  homogeneous  layers.  It  is  able  to  include  all  phases 
resulting  from  the  interaction  of  the  propagating  wavefield  with  all 

interfaces  in  the  model  and  the  free  surface.  It  is  thus  ideal  as  a  tool  to  predict 

the  likely  features  of  coda  generated  simply  by  layer  resonance.  We  have 
used  the  constant  Q  Ngendei  model  proposed  by  Sereno  (1986).  This  model  is 

distinguished  by  a  5.5  km  thick  water  column  and  a  low  velocity  (a=1.6  km/s, 
P=0.116  km/s)  70  m  thick  sediment  layer. 

SYNTHETIC  RESULTS 

We  have  computed  synthetic  Green's  functions  corresponding  to  a  lateral 
offset  between  the  source  and  receivers  of  1000  km.  To  carry  out  the  actual 
computation  two  separate  sources,  each  consisting  of  vertical  and  horizontal 

point  forces,  were  located  at  the  sediment/water  interface  and  0.5  m  above  the 
interface  in  the  water  column.  The  receiver  was  located  at  15  km  depth.  Using 
reciprocity  we  obtain  the  response  to  a  source  buried  at  15  km  when  observed 
at  the  shallower  depths  (Shaw  and  Orcutt,  1984).  These  receiver  locations 
enable  us  to  predict  the  relative  abilities  of  hydrophones  and  seafloor 
seismometers  for  sensing  the  resonance  of  energy  in  the  sediment  layer  and 
water  column. 

In  figure  1  we  display  the  vertical  normal  stresses  due  to  a  vertical  point  force 
(£zz0)  at  both  receiver  depths.  We  expect  this  Green's  function  should  be 
strongly  influenced  by  the  propagation  of  high  phase  velocity  comprcssional 
energy.  Frequency  and  time  representations  of  these  functions  arc  displayed 
in  the  upper  and  lower  halves  of  the  figure  respectively.  The  portions  of  the 
time  series  used  in  computing  the  spectra  arc  indicated  by  horizontal  bars  in 
the  lower  figures.  The  spectral  estimates,  computed  using  a  multi-taper 
algorithm  (Thompson,  1982),  consist  of  the  sum  of  5  sub-estimates,  computed 


using  a  time-bandwidth  product  of  3.0.  The  upper  and  lower  spectra, 

corresponding  to  the  receivers  at  5.5495  and  5.5000  km  depth  respectively, 
both  are  dominated  by  a  regular  modulation  with  a  spacing  of  0.136  Hz.  This 
modulation  was  explained  by  Sereno  (1986)  as  being  due  to  the  resonance  of 
acoustic  energy  in  the  water  column.  The  predicted  eigenfrequencies 
corresponding  to  this  resonance  are  indicated  by  the  vertical  dashed  lines. 
There  is  a  striking  correlation  between  synthetically  computed  and  predicted 
spectral  peaks.  It  is  clear  that  this  Green's  function  is  strongly  influenced  by 
the  resonance  of  acoustic  energy  in  the  water  column.  By  examining  the  time 
series  there  is  clearly  a  periodicity  of  roughly  7  seconds  (the  two  way  travel 

time  of  acoustic  energy  in  the  water  column)  following  the  onset  of  Pn  energy 
about  125  seconds  after  the  event.  From  this  exercise  we  conclude  that  sensors 
on  the  seafloor,  and  above  in  the  water  column,  should  be  sensitive  to  the 

resonance  of  acoustic  energy  in  the  water  column. 

In  figure  2  we  display  the  same  Green's  function,  however  the  time-window 
considered  was  shorter  (20  sec)  and  restricted  to  the  time  immediately 

following  the  Pn  onset.  We  expect  that  this  Green's  function  should  also  be 
influenced  by  compressional  energy  resonance  in  the  sediment  layer.  The 

spectra  and  time  series  in  this  figure  are  arranged  in  the  same  manner  as 
those  in  figure  1,  however  the  spectral  estimates  were  made  using  7  sub¬ 

estimates  computed  using  a  time-bandwidth  product  of  4.0.  We  observe  a  broad 
power  high  centered  close  to  the  eigenfrequency  predicted  for  the  resonance 
of  compressional  energy  in  the  sediment  layer.  The  spectral  rolloff  occurring 
above  8  Hz  is  due  to  the  taper  applied  when  converting  the  green's  functions 

from  the  frequency  domain  to  the  time  domain  and  is  not  due  to  resonance. 
The  predicted  frequency  spacing  of  this  modulation  is  11.4  Hz  and  thus  the 

synthetics  computed  to  a  nyquist  frequency  of  only  9  Hz  only  illustrates  the 

first,  broad,  spectral  peak. 

In  figure  3  we  illustrate  the  same  results  for  the  vertical  normal  stress  due  to  a 
horizontal  point  force  (£zzi).  We  expect  this  Green's  function  should  be 

greatly  influenced  by  radial  particle  motion  and  thus  should  be  strongly 

dependent  on  the  resonance  of  shear  energy  in  the  sediment  layer.  We  have 
indicated  the  frequencies  of  the  predicted  spectral  peaks  by  the  dashed 

vertical  lines  in  the  upper  half  of  the  figure  (the  frequency  spacing  is  0.83 
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Hz).  We  have  computed  spectra  from  90  seconds  of  the  Sn  coda.  The  lower 
spectrum  was  computed  from  the  receiver  located  at  the  sediment/water 
interface.  As  expected  it  is  strongly  influenced  by  shear  wave,  sediment  layer, 
resonance.  The  green's  function  in  the  water  column,  just  0.5  m  above,  is 
clearly  oblivious  to  this  resonance.  The  spectrum  is  influenced  by  acoustic, 
water  layer,  resonance  and  little  else.  By  inspection  of  the  Green's  functions 
in  the  time-domain,  the  water  multiples  are  obvious  in  the  water  column  but 
much  inter-multiple  energy  due  to  sediment  layer  resonance,  that  is  present 
in  the  lower  time  series,  is  absent.  We  believe  there  is  a  simple  explanation  for 
this.  In  figure  4  we  display  the  modulus  of  the  transmission  coefficient  for 
upward  traveling  energy,  converting  from  shear  to  compressional 
propagation  at  the  sediment/water  interface.  Although  the  transmission 
coefficient  is  a  monotonically  increasing  function  of  ray  parameter,  it 
remains  below  0.007  over  the  range  of  slownesses  we  expect  to  dominate  the  Pn 
and  Sn  coda  (0.12  to  0.35  s/km).  As  one  would  intuitively  expect,  considering 
the  particle  motions  involved,  shear  energy  does  not  readily  convert  to 
compressional  energy  at  the  boundary  between  these  low  velocity  horizons. 
We  conclude  that  the  WIH  dataset  is  inappropriate  when  seeking  evidence  of 
shear  wave  resonance  in  the  sediment  layer  and  we  do  not  consider  it  further. 

RECORDED  RESULTS 

We  now  consider  actual  coda  recorded  by  the  Wake  Island  Hydrophones  and 
the  Ocean  Bottom  Seismometers  and  Hydrophones  deployed  during  the  Ngendei 
expedition.  In  figure  5  we  display  the  spectra  and  time  series  corresponding  to 
an  event  roughly  1000  km  from  the  receivers.  We  indicate  by  the  dashed 
vertical  lines  the  predicted  eigenfrequencies  corresponding  to  compressional 
energy  resonating  in  the  sediment  layer  (assuming  a  net  it  phase  shift 
between  successive  multiples).  The  lower  spectrum,  computed  from  50  seconds 
of  Pn  coda  recorded  by  the  vertical  component  seismometer  (pictured  at  the 
bottom),  displays  a  clear  modulation  consistent  with  the  prediction  at 
frequencies  below  40  Hz  (the  point  at  which  signal  falls  into  noise).  The 
spectrum  computed  from  the  same  portion  of  the  Pn  coda  recorded  by  the 
hydrophone  displays  a  similar,  though  less  evident,  modulation;  however,  the 
positions  of  the  power  highs  suggest  that  successive  compressional  sediment 
layer  multiples  are  in  phase.  The  power  highs  are  consistent  with  the  second, 
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solid,  set  of  vertical  lines  computed  assuming  no  net  phase  shift.  This 
phenomenon  is  not  restricted  to  this  event  but  is  common  to  all  we  have 
examined  that  have  induced  the  resonance  of  compressional  energy  in  the 
sediment  layer.  It  is  possible  to  induce  compressional  energy  resonance  in  the 
sediment  layer  by  introducing  energy  from  above,  or  below.  As  is  discussed  by 
Sereno  (1986),  arrivals  from  below  are  followed  by  an  infinite  sequence  of 
echos  each  n  phase  shifted  from  its  immediate  neighbors  in  time.  The 
cumulative  modulus  of  reflection  incurred  by  bouncing  at  the  sediment/water 
and  then  the  sediment/crust  interfaces  is,  assuming  the  constant  Q  Ngendei 
model,  0.029.  The  primary  arrival  (with  an  amplitude  of  1.0)  is  followed  by  a 
second  with  an  amplitude  of  0.029,  and  a  third  of  .00084  and  so  on  -  the  echos 
die  away  quickly.  Significant  resonances  should  be  induced  by  downward 
traveling  energy  that  has  reflected  at  the  free  surface.  The  primary 
downward  traveling  arrival  will  pass  the  sediment  water  interface  and  be 
followed  by  an  echo  after  the  energy  reflects  at  the  sediment/crust  interface. 
Energy  propagating  with  a  phase  velocity  of  8.2  km/s  undergoes  no  phase 
shift  either  in  transmitting  into  the  sediment  layer,  or  by  the  reflection  at  the 
top  of  the  crust,  and  thus  the  primary  arrival  and  its  echo  should  be  in  phase. 
The  resultant  spectral  modulation  will  have  peaks  at  multiples  of  the 
frequency  of  the  first  (see  the  solid  vertical  lines  in  figure  5).  The  echo 
should  have  an  amplitude  of  0.79  relative  to  the  initial  arrival.  This  form  of 
resonance  is  nearly  truncated  after  the  first  echo.  All  subsequent  echos  have 
insignificant  amplitudes  relative  to  the  first.  For  this  theory  to  be  tenable  the 
observed  "anti-correlation"  between  the  hydrophone  and  seismometer  spectra 
must  not  be  present  in  the  portion  of  the  coda  after  the  Pn  onset  and  before 
the  arrival  of  the  first  multiple  off  the  free  surface  (a  span  of  7.3  seconds).  In 
figure  6  we  look  at  the  vertical  component  seismometer  and  hydrophone 
spectra  computed  from  the  initial  6.5  seconds  of  coda  after  the  Pn  onset  of  an 
event  recorded  by  OBS  karen.  Virtually  no  broad  modulation  can  be  discerned 
in  either  of  these  spectra  (there  is  only  a  slight  suggestion  of  spectral  peaks  at 
5.7  and  17.1  Hz),  however  there  is  a  clear  coincidence  of  smaller  scale  spectral 
features.  When  we  look  at  the  time  span  from  7.0  to  13.5  seconds  after  the  Pn 
onset  (figure  7)  we  see  that  the  modulation  is  much  more  apparent  and  a 
distinct  anti-correlation  between  seismometer  and  hydrophone  spectra  is 
present.  In  this  time  span  the  energy  from  the  free  surface  has  had  enough 
time  to  return  to  the  seafloor.  There  is  likely  a  competition  between  two  types 
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of  compressional  sediment  layer  resonances.  One  that  includes  an  infinite 
sequence  of  quickly  decaying  echos  that  are  n  phase  shifted  from  their 

neighbors,  and  the  other  consisting  of  only  a  single  significant  echo  that  is  in 
phase  with  the  primary  arrival.  It  appears  that  in  the  Ngendei  dataset  the 

latter  form  of  resonance  is  dominant  at  the  level  of  the  hydrophone  and  the 

former  dominates  the  seismometer  -  coupled  to  the  sediment  layer.  In  the 

synthetics  it  appears  that  the  former  type  of  resonance  dominates  both  the 

sensor  in  the  water  column  and  the  one  at  the  interface. 

The  event  in  figure  8  was  recorded  by  the  WIH  array.  It  was  an  event  with  mb 
of  6.0  and  occurred  at  a  depth  of  159  km  and  a  range  of  2170  km.  The  spectra 
display  signal  and  noise  levels  and  were  computed  from  the  portions  of  the 

time  series  delineated  by  the  horizontal  bars.  The  time  series  obtained  from 
each  of  the  6  seafloor  sensors  that  recorded  this  event  were  stacked  while 

beamfoiming  for  the  Pn  phase.  Signal  remains  above  noise  in  the  Pn  coda 
produced  by  this  event  to  a  frequency  of  22  to  25  Hz.  A  subtle,  broad, 

modulation  is  evident  at  lower  frequencies  and  corresponds  to  the 
eigenfrequencies  predicted  assuming  no  net  phase  shift  between  successive 
water  multiples.  The  spacing  of  these  power  highs  suggests  a  two  way  travel 
time  in  the  sediment  layer  of  roughly  0.15  seconds. 

In  figure  9  we  display  the  same  event  pictured  in  figure  5.  We've  magnified 
the  frequency  scale  of  the  upper  half  of  the  figure  to  illustrate  the  fine  scale 
modulation  consistent  with  resonance  of  acoustic  energy  in  the  water  column 
(see  the  vertical  dashed  lines).  Both  spectra  have  a  fine  modulation;  however, 
the  upper  spectrum,  computed  from  the  hydrophone  channel,  possesses  a 
more  pronounced  modulation  at  this  scale.  We  see  a  similar  result  when 
examining  a  second  event  recorded  by  the  WIH  array  (figure  10).  This  event 
was  an  mb  5.4  earthquake,  33  km  deep  which  occurred  at  a  range  of  2057  km 
from  the  WIH  array. 


CONCLUSIONS 

We  have  observed  significant  spectral  modulation  present  in  the  coda  of 
events  recorded  by  seafloor  seismometers  and  hydrophones  and  near-seafloor 
sensors  in  the  Wake  Island  Hydrophone  array.  These  modulations  are 
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pronounced  and  consistent  with  those  present  in  synthetic  spectra  in  which 
we  know  the  regular  spectral  peaks  are  due  to  layer  resonance  only.  Our 
results  substantiate  the  conclusions  made  by  Sereno  and  Orcutt  (1985a,b)  that 
layer  resonance  is  a  significant,  perhaps  dominant,  source  of  coda  in  the 
oceanic  environment. 
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FIGURE  CAPTIONS 


Figure  1.  Synthetic  Green's  functions  (IZz0)  and  corresponding  frequency 
power  spectra  at  depths  of  5.5495  and  5.5000  km  due  to  a  vertical  point  force 
buried  15  km  at  a  range  of  1000  km.  In  figures  1,2,3,5,6,7,9  and  10  the  upper 
spectral  estimates  correspond  to  the  upper  time  series.  In  each  case  the 
estimates  were  computed  from  the  portions  of  the  seismograms  indicated  by 
the  horizontal  bars.  The  spectral  estimates  were  computed  using  a  time- 
bandwidth  product  of  3. 

Figure  2.  Synthetic  Green's  functions  (EzzO)  and  corresponding  frequency 

power  spectra  at  depths  of  5.5495  and  5.5000  km  due  to  a  vertical  point  force 
buried  15  km  at  a  range  of  1000  km.  The  spectral  estimates  were  computed 

using  a  time-bandwidth  product  of  4. 

Figure  3.  Synthetic  Green's  functions  (ZZzl)  and  corresponding  frequency 

power  spectra  at  depths  of  5.5495  and  5.5000  km  due  to  a  horizontal  point  force 
buried  15  km  at  a  range  of  1000  km.  The  spectral  estimates  were  computed 

using  a  time-bandwidth  product  of  3. 

Figure  4.  The  modulus  of  the  transmission  coefficient  corresponding  to  the 
conversion,  at  the  sediment/water  interface,  of  shear  energy  propagating 
upward  and  converting  to  upward  propagating  compressional  energy. 

Figure  5.  Power  spectral  estimates  and  seismograms  corresponding  to  an 
earthquake  1000  km  from  the  sensor.  The  data  were  collected  by  an  Ocean 
Bottom  Seismometer  and  Hydrophone  deployed  in  the  S.W.  Pacific.  The  vertical 
dashed  and  solid  lines  indicate  the  predicted  locations  of  spectral  peaks  due  to 
compressional  energy  resonance  in  the  sediment  layer  where  successive 

multiples  are  phase  shifted  by  n  and  0.0  respectively. 

Figure  6.  In  this  figure  we  consider  the  event  displayed  in  figure  5; 

however,  we  only  consider  the  first  6.5  seconds  of  coda,  following  the  Pn  onset, 
in  the  spectral  estimates.  The  solid  vertical  lines  indicate  the 
eigenfrequencies  associated  with  the  resonance  of  compressional  energy  in 

the  sediment  layer  in  which  successive  multiples  are  n  phase  shifted. 

Figure  7.  In  this  figure  we  consider  the  event  displayed  in  figure  5; 

however,  we  only  consider  the  portion  of  the  coda  arriving  from  7.0  to  13.5 

seconds  after  the  Pn  phase.  The  dashed  and  soiid  vertical  lines  indicate  the 

eigenfrequencies  associated  with  the  resonance  of  compressional  energy  in 

the  sediment  layer  in  which  successive  multiples  are  phase  shifted  by  n  and 
0.0  respectively. 

Figure  8.  Signal  and  noise  power  spectral  estimates  and  seismogram 

corresponding  to  an  mb  6.0  earthquake  which  occurred  at  a  depth  of  159  km 
and  a  range  of  2168  km  from  the  Wake  Island  Hydrophone  array.  The  vertical 
dashed  lines  indicate  the  predicted  locations  of  spectral  peaks  due  to 

compressional  energy  resonance  in  the  sediment  layer  where  successive 

multiples  are  in  phase. 
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Figure  9.  Power  spectral  estimates  and  seismograms  corresponding  to  the 

event  displayed  in  figure  5.  The  vertical  dashed  lines  indicate  the  predicted 
locations  of  spectral  peaks  due  to  compressional  energy  resonance  in  the 

water  column. 

Figure  10.  Power  spectral  estimate  and  seismogram  corresponding  to  mb  5.4 
earthquake  which  occurred  at  a  depth  of  33  km  and  a  range  of  2057  km  from 

the  Wake  Island  Hydrophone  array.  The  vertical  dashed  lines  indicate  the 

predicted  locations  of  spectral  peaks  due  to  compressional  energy  resonance 
in  the  water  column. 
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REGIONAL  CONTINENTAL  PROPAGATION 


Michael  A.H.  Hedlin  and  John  A.  Orcutt 


103 


SUMMARY 


We  are  currently  examining  data  collected  in  central  Asia  and  Scandinavia  to 

gain  a  deeper  insight  into  the  nature  of  the  propagation  of  seismic  phases 

through  the  continental  crust.  As  has  been  discussed  in  the  previous  chapter, 
two  popular  theories  exist  that  attempt  to  explain  the  mechanism  by  which 
seismic  coda  is  generated  in  the  oceanic  crust.  These  theories  involve  the 
scattering  of  seismic  energy  by  small  scale  inhomogeneities  (Richards  and 
Menke,  1983;  Menke  and  Chen,  1984)  and  the  resonance  of  energy  in  low 
velocity  horizons  (Sereno  and  Orcutt,  1985a, b).  We  attempt  to  determine  the 
likely  dominant  coda  generation  mechanism  present  in  the  continental  crust 
by  processing  the  data  in  two  ways.  We  compute  frequency  versus  time 

expansions  of  the  energy  in  synthetic  and  actual  coda  in  a  manner  described 
by  Hedlin,  Minster  and  Orcutt  (1989).  We  also  compute  frequency  versus 
wavenumber  spectra  from  isolated  phases  in  the  synthetic  and  real  coda.  The 
synthetic  seismograms  are  computed  using  Wavenumber  Integration  (Luco 

and  Apsel,  1983;  Apsel  and  Luco,  1983). 
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FREQUENCY-TIME  ANALYSIS 


We  have  found  frequency-time  display  of  seismograms,  known  as  sonograms 
in  acoustics  (Market  and  Gray,  1976),  to  be  particularly  useful  for  coda  analysis 
because  the  time  dependence  information  is  preserved.  We  have  first  applied 
this  processing  technique  to  a  synthetic  data  set  created  using  Wavenumber 
Integration.  This  algorithm  requires  that  the  Earth  be  modeled  as  consisting 
of  a  stack  of  homogeneous  layers  and  includes  all  seismic  phases  resulting 
from  the  interaction  of  the  propagating  wavefield  with  the  interfaces  within 
the  model  and  the  free  surface.  It  is  thus  ideal  as  a  tool  to  predict  the  likely 
features  of  coda  generated  in  the  absence  of  scattering,  simply  by  the 
resonance  of  energy  in  low  velocity  horizons.  As  was  demonstrated  by  Hedlin 
et  al  (1988)  it  is  possible  to  generate  significant  coda  u~ing  an  Earth  model 
consisting  of  homogeneous  layers  and  no  water  column.  The  model  used  was 
representative  of  the  continental  crust  in  the  vicinity  of  the  Semipalatinsk 
nuclear  test  site  in  central  Asia  and  was  capped  by  a  low  velocity  horizon 
representing  a  near-surface  weathered  layer  or  sediment.  In  figures  la  and 
lb  are  displayed  synthetic  seismograms,  and  accompanying  sonograms,  of  the 
vertical  and  radial  components  of  motion  caused  by  an  explosion  buried  30  m 
at  a  range  of  400  km.  Since  the  Earth  model  consists  of  homogeneous  layers, 
the  synthetic  coda  consists  solely  of  reverberations  primarily  concentrated  in 
the  surface  layer.  Sereno  and  Orcutt  (1985a)  showed  that  coda  produced  by 
layer  resonance  is  enriched  in  power  at  certain  equispaced  eigenfrequencies. 
Peak  spacing  is  inversely  proportional  to  the  two-way  vertical  travel  time 
across  the  layer,  and  peak  width  is  controlled  by  the  impedance  contrasts 
between  the  layer  and  its  neighbors  as  well  as  the  attenuative  properties  of 
the  layer  itself.  As  was  discussed  by  Hedlin  et  al  (1988),  and  is  apparent  in  this 
figure,  when  the  receiver  is  a  significant  distance  from  the  source  (so  that  the 
incident  energy  has  high  phase  velocity)  the  period  of  the  modulation 
observed  on  vertical  or  horizontal  component  sensors  is  largely  independent 
of  time  in  the  coda.  The  eigenfrequency  spacing  in  figure  lb  is  closer  than  in 
figure  la,  because  at  this  range  the  radial  component  i«  dominated  by  shear 
energy,  and  the  vertical  component  by  faster  compressional  energy. 
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Wc  expect  that  actual  coda  produced  largely  by  the  resonance  of  energy  in 
low-velocity  strata  should  have  frequency-time  structure  similar  to  that 
possessed  by  these  synthetics.  We  have  examined  two  independent  data  sets  in 
search  of  these  time  independent  spectral  features.  The  first  data-set  consists 
of  high  frequency  (250  samples/s)  recordings  of  regional  mine  explosions  and 
single-event  calibration  explosions  made  by  a  triangular  array  of 
seismometers  deployed  in  Kazakhstan  U.S.S.R.  (Eissler  et  al,  1988;  Thurber  et  al, 
1988).  The  second  data-set  consists  of  recordings  of  regional  mine  explosions 
and  earthquakes  digitized  at  40  samples/s  by  the  NORESS  array  in  Norway 
(Sereno  et  al,  1987).  As  has  been  discussed  by  Hedlin,  Minster  and  Orcutt  (1989) 
and  Hedlin  et  al  (1989)  time  independent  spectral  features  have  been  observed 
in  a  number  of  events  in  central  Asia  and  Scandinavia  but  these  events  have 
all  been  identified  as  being  mine  explosions.  The  time-independent  features 
are  believed  to  originate  at  the  source  from  a  detonation  practice  known  as 
ripple  firing,  in  which  a  number  of  sub-explosions  are  closely  grouped  in 
space  and  time.  In  figure  2  we  display  the  sonogram  corresponding  to  a  Blasjo 
mine  explosion  recorded  at  a  range  of  302  km  by  the  NORESS  array.  All  events 
believed  to  be  either  calibration  explosions  or  earthquakes  (and  thus  likely  to 
have  relatively  uncolored  source  spectra)  give  rise  to  coda  apparently  free  of 
time-independent  spectral  features.  In  figures  3  and  4  we  display  the 
sonograms  corresponding  to  two  earthquakes  which  occurred  in  Scandinavia 
and  were  recorded  by  the  NORESS  array.  We  conclude  that  the  continental 
crust  is  not  giving  rise  to  significant  layer  resonance  and  thus  the  dominant 
coda  generation  mechanism  in  the  vicinity  of  these  two  sets  of  seismometers 
likely  involves  some  other  process.  The  most  likely  alternative  is  the 
scattering  of  energy  by  inhomogeneities.  The  minor  time-independent 
feature  present  in  the  coda  of  the  event  displayed  in  figure  4  exists  prior  to  the 
onset  and  is  clearly  a  long  lived  spectral  line  in  the  noise  unrelated  to  the 
event. 

FREQUENCY-WAVENUMBER  ANALYSIS 

We  feel  we  can  gain  additional  insight  regarding  the  source  of  seismic  coda  in 
the  continental  crust  by  isolating  various  phases  within  seismic  coda  and 
exploiting  array  data  to  compute  frequency  versus  wavenumber  spectra.  It 
has  been  demonstrated  by  many  authors  (including  Mykkeltveit  and  Bungum, 
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1984;  Dainty,  198S)  that  wavenumber  spectra  allow  the  simultaneous  estimation 
of  phase  arrival  azimuth  and  velocity.  For  a  preliminary  synthetic  result  we 
have  used  seismograms  computed  using  the  Earth  model  discussed  in  the 
previous  section  and  linearly  stacked  it  upon  itself  after  time  shifting  for  a 
phase  coming  from  a  back  azimuth  of  240°  with  a  phase  velocity  of  8.2  km/s 
assuming  sensor  locations  in  the  NORESS  array  (figure  S).  In  this  example  the 
source  was  placed  300  km  from  the  array.  In  figure  6  we  display  the 
wavenumber  spectrum  computed  from  this  synthetic  event  considering  only 
the  energy  immediately  after  the  compressional  onset  in  the  frequency  band 
from  4.5  to  5.0  Hz.  This  synthetic  event  defines  the  extreme  limit  in  resolution 
of  the  NORESS  array.  Because  of  the  manner  in  which  we  constructed  the 
input,  all  phases  arrive  with  the  same  phase  velocity  and  from  the  same 
azimuth.  As  a  result  the  true  wavenumber  spectrum  is  a  delta  function  at  a 
back  azimuth  of  240°  and  a  Ikl  of  0.58  cycles/km.  In  figure  6  it  is  clear  that  the 
center  of  the  energy  distribution  is  at  the  correct  location;  however,  the 
impulse  response  of  the  array  has  smeared  the  actual  wavenumber  field  out  to 
a  "patch"  with  a  half-width  of  roughly  0.4  cycles/km.  In  the  next  section  we 
will  discuss  the  implications  of  the  impulse  response  somewhat  further.  We 
expect  that  the  wavenumber  spectrum  of  a  wavefield  that  has  propagated 
through  a  medium  solely  in  the  plane  joining  the  source  and  receiver  should 
resemble  this  event  to  the  extent  that  the  energy  should  be  confined  to  this 
plane.  Wavenumber  spectra  computed  from  the  coda  of  a  Norwegian  quarry 
blast  (at  a  range  of  302  km,  and  a  back  azimuth  of  240°  from  NORESS)  is  shown 
in  figures  7  to  10.  In  figure  7  we  display  the  wavenumber  spectrum  computed 
from  energy  in  the  frequency  band  from  1.5  to  5.0  Hz  arriving  in  the  window 
from  the  compressional  onset  to  6.5  seconds  later.  At  this  range  we  expect  this 
window  to  include  crustal  phases  as  well  as  the  Pn  phase.  The  energy  is  fairly 
well  constrained  to  lie  on  the  plane  joining  the  source  and  the  array.  In 
figure  8  we  include  an  additional  13.5  seconds  of  coda.  There  is  now  a  fairly 
considerable  amount  of  apparent  "off  axis"  energy.  We  have  not  resolved  the 
ambiguity  over  whether  this  indicates  out  of  plane  propagation  or  is  due  to 

coherent  noise.  The  average  signal  levels  in  this  window  are  somewhat  lower 
than  those  present  in  the  previous  one.  In  figures  9  and  10  we  consider  the  Sn 

onset  and  the  following  Lg  coda.  Again  there  is  considerable  energy  evident 

due  to  coherent  noise  or  out  of  plane  propagation.  We  are  currently 

examining  additional  events,  using  a  sliding  window  wavenumber  analysis,  in 
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an  attempt  to  establish  the  origin  of  the  apparent  "off-axis"  energy.  In 
addition  we  are  producing  a  full  suite  of  synthetic  Green's  functions  for  the 
NORESS  array  by  considering  source  to  receiver  offsets  consistent  with  each 
individual  element  in  the  array. 


We  can  estimate  the  ability  of  an  array  to  resolve  features  in  the  wavenumber 
domain  by  computing  the  impulse  response  function,  or  the  response  of  the 
array  to  the  simultaneous  excitation  of  every  sensor  with  a  delta  function.  We 
define  the  input  to  the  array  as: 


N 


u  (x,y)  =  5^  5(x-x.)5(y-y.) 

j=l 


(1) 


where  N  is  the  number  of  sensors  and  (xj.yj)  are  the  locations  of  the  sensors. 
The  response  in  the  wavenumber  domain  is  given  by 

XN  N 


U(kx,ky)  =  £XXX  8(x-x.)S(y-y.)e 

x=xjy=yj  j=i 


i2rt(kxx+kvy) 


(2) 


where  (kx,ky)  and  the  x  and  y  wavenumbers  in  cycles/km.  The  triple  sum 
simplifies  to: 


1  ^  -i2n(k  x.+k  y.) 

U(k«V  =  N  X  e 

j=l 

The  amplitude  of  the  impulse  response  equals: 


itxk>i.&  /  J ;i2’"k,t');£  ei2”<k-r") 

V  j=l  m-1 


(3) 


(4) 


where  k  =  (kx,ky).  This  can  be  simplified  to: 

I  U(k) 

The  impulse  response  of  the  NORESS  array  is  shown  in  figure  11. 


V  m=l  j=l 


(5) 

When  spatial 
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periodogram  spectral  estimates  are  computed  (ie  no  spatial  taper  is  applied  to 
the  recorded  wavefield)  the  true  wavenumber  spectrum  of  an  event  is 
convolved  with  this  function.  As  we've  discussed  earlier,  the  half-width  of  the 
main  lobe  (roughly  the  aperture  of  the  array)  is  approximately  0.4  cycles/km 
and  defines  the  limit  of  resolution  in  the  wavenumber  domain.  No  two  events 
separated  by  less  than  the  resolution  can  be  resolved  as  being  distinct  events. 

Further  corruption  of  the  computed  wavenumber  field  occurs  because  of  the 

sidelobes  which,  for  this  array,  are  small  but  not  infinitesimally  so.  The  true 
wavenumber  spectrum  has  been  blurred  to  a  circular  "patch"  with  a  half 

width  of  0.4  cycles/km.  The  NORESS  array  is  extremely  well  designed  but  we 
feel  that  some  improvement  in  resolution  can  be  gained  by  the  judicious 
application  of  spatial  tapers.  Currently  we  are  extending  the  theory  of 
Thompson  (1982),  who  dealt  with  the  computation  of  optimal  tapers  for  evenly 
sampled  1-D  time  series,  to  the  2-D  array  data.  The  primary  obstacle  to  the 
application  of  the  method  is  the  unevenness  of  the  spatial  sampling  which 
complicates  the  computation  of  the  eigentapers.  The  theory  of  Thompson 
(1982)  is  readily  applicable  to  the  time-frequency  transformation  required  by 

this  method.  As  is  discussed  by  Hedlin,  Minster  and  Orcutt  (1989)  it  will  help  us 
isolate  particular  phases  by  allowing  shorter  time  windows. 

CONCLUSIONS 

Our  preliminary  analysis  of  the  evolution  of  frequency  spectra  as  a  function 
of  time  in  seismic  coda  suggests  that  layer  resonance  is  not  playing  a  major 
role  in  the  generation  of  seismic  coda  in  the  continental  crust.  We  are 

currently  employing  an  analysis  of  the  evolution  of  the  wavenumber 

spectrum  as  function  of  time  in  an  attempt  to  determine  if  the  likely 

alternative  mechanism,  the  scattering  of  energy  by  small  scale 

inhomogeneities,  is  responsible. 
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FIGURE  CAPTIONS 


Figure  1.  Synthetic  vertical  and  radial  acceleration  sonograms  and 
corresponding  time  series  corresponding  to  an  explosive  source  buried  30  m  at 
a  range  of  400  km.  The  model  used  in  this  calculation  consisted  of  a  stack  of 
homogeneous  layers  above  a  half  space.  The  model  is  capped  by  a  100  m  thick 
low  velocity  horizon  representing  a  weathered  or  sediment  layer. 

Figure  2.  Seismogram,  and  binary  sonogram,  corresponding  to  a  Blasjo 

mine  explosion  recorded  at  a  range  of  302  km  by  the  NORESS  array. 

Figure  3.  Seismogram,  and  binary  sonogram,  corresponding  to  a  presumed 

earthquake  recorded  at  a  range  of  343  km  by  the  NORESS  array. 

Figure  4.  Seismogram,  and  binary  sonogram,  corresponding  to  a  presumed 

earthquake  recorded  at  a  range  of  513  km  by  the  NORESS  array. 

Figure  5.  The  NORESS  array. 

Figure  6.  The  wavenumber  spectrum  of  a  synthetic  event  at  a  range  of  300 
km  produced  using  the  model  described  in  the  caption  for  figure  1.  The 
recording  by  the  NORESS  array  was  simulated  by  placing  the  displayed  time 
series  at  the  locations  of  the  sensors  in  the  array  and  offsetting  in  time  for  a 
phase  arriving  from  an  azimuth  of  240°  with  an  apparent  velocity  of  8.2  km/s. 

Figure  7.  The  wavenumber  spectrum  of  a  Blasjo  mine  explosion.  Only  the 

portion  of  the  Pn  coda  indicated  by  the  lower  figure  was  considered.  The 

estimate  was  constructed  from  sub-estimates  spanning  the  frequency  range 

from  1.5  to  5.0  Hz.  The  time  series  has  been  beamformed  for  the  Pn  phase. 

Figure  8.  The  wavenumber  spectrum  of  a  Blasjo  mine  explosion.  Only  the 

portion  of  the  Pn  coda  indicated  by  the  lower  figure  was  considered.  The 

estimate  was  constructed  from  sub-estimates  spanning  the  frequency  range 

from  1.5  to  5.0  Hz.  The  time  series  has  been  beamformed  for  the  Pg  phase. 

Figure  9.  The  wavenumber  spectrum  of  a  Blasjo  mine  explosion.  Only  the 

portion  of  the  Sn  coda  indicated  by  the  lower  figure  was  considered.  The 

estimate  was  constructed  from  sub-estimates  spanning  the  frequency  range 

from  1.5  to  5.0  Hz.  The  time  series  has  been  beamformed  for  the  Sn  phase. 

Figure  10.  The  wavenumber  spectrum  of  a  Blasjo  mine  explosion.  Only  the 

portion  of  the  Sn  coda  indicated  by  the  lower  figure  was  considered.  The 

estimated  was  constructed  from  sub-estimates  spanning  the  frequency  range 
from  1.5  to  5.0  Hz.  The  time  series  has  been  beamformed  for  the  Lg  phase. 

Figure  11.  The  impulse  response  of  the  NORESS  array. 


Ill 


Figures  la  and  lb.  Synthetic  vertical  and  radial  acceleration  sonograms.  The  source  is  an 
explosion  buried  30  meters  at  a  range  of  400  km. 
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